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SYMDPSIS 

CAIMOMCAL STABILITY MA.TRICES FOE LINEAR TPIE- 
invariant systei/b and THEIE AEPLICATIOMS 

Ramamoortiiy Viswanatiian 

Sciiwarz in 1956 introduced a method of checking 
the stability of a linear time-invariant system 
described by the vector differential equation 

X = A X 

This method requires transforming the matrix a through 
a series of elementary transformations to a canonical 
matrix B given by 



continued fractions, Schwarz 


Using the concept of 



SYNOPSIS 


CANONICAL STABILITY miRICES FOE LINEAR TB!E- 
INYARIANT SYSTEI/B Al® THEIR APPLICATIONS 

Ramartioortiiy Viswanatliaa 


Schwarz in 1956 Introduced a method of checking 
the stability of a linear time-invariant system 
described by the vector differential equation 

z = A X 


This method requires transforming the matrix a through 
a series of elementary transformations to a canonical 
matrix B given by 


0 10 
-b, 0 1 

0 - 1 .^ 0 


0 0 
0 0 
0 0 


0 

0 

Using the concept 


0 0 

0 0 

of continued 


0 1 



— r 


fractions, Schwarz 



xi 


showed -thsrt; the given system, is asyj^jtotically stable if and 
only if all the h^’ s , i=l ,S , . . . ,n of B are positive 
and that the number of stable eigenvalues is equal to 
the number of positive terms in the sequence b^^, 

Recently in the control 
systems literature there has been a great deal of 
interest shown towards finding uses for the above 
Schwarz canonical matrix. Through this canonical matrix, 
Liapunov function for the system can be quite simply 
constructed and a certain class of performance measures 
can be evaluated. It has been used to establish a link 
between the second method of Liapunov and the Routh- 
Hurwitz stability critorlon. 

C In this thesis, three different- but 
forms of the Schwarz matrix are identified and the 
similarity transformations between them are determined 
The problem of transforming a given system described by i 

X = X 

to one in Schwarz form 

y = B y . 

I 

is dealt with systematically. When A is in phase-' 
variable canonical form, the associated transformation'! 
is shown to be quite simple and it requires the; 



xii 


knowledge of tiie elements of tke Routii-array of the 
system. The transformation between the Schwarz form 
and its Jordan form is also determined and it turns out 
to be quite analogous to the one between phase-variable 
form and its Jordan form.’ Schwarz canonical form for 
controllable multivariable systems Is proposed i- This 
consists of a set of coupled subsystems each of which 
being in Schwarz form, 'transforming any arbitrary 
controllable multivariable system to this canonical 
form is also discussed. 

A network interpretation is given to the three 
forms of the Schwarz matrix. By a proper choice of 
state variables, the state equations of a resistlvely 
terminated L-C ladder network are shown to be in any 
one of these forms. Also, given a system in Schwarz 
form, a procedure is given to evaluate the components 
of the associated resistlvely terminated L-0 ladder 
network. Thus such a network interpretation is shown 
to be useful both for analysis and synthesis problems. 
Network interpretation for multivariable systems in 
Schwarz canonical form is also dealt with. 

An important use of the Schwarz canonical form 
in simplifying large dynamic systems is also investiga- 
ted in this thesis. Problems involving large dynamic 
systems demand complex computational schemes and thus 



xiii 


it Is desirable to make reasonable approximations to 
simplify such large systems. The simplification method 
proposed in this thesis first transforms the given 
’large’ system to one in Schwarz form, compares the 
elements bj_’s of the Schwarz form and thus determines 
whether any simplification is possible or not and 
finally carries out the actual simplification. Several 
examples are considered for illustrating the method. 
Some difficulties encountered by this method are also 
pointed out. The major advantage of this method is that 
it does not require the computation of the system eigenr 
values and eigenvectors. 

The Schwarz matrix is used to derive a new cano- 
nical matrix for linear time- invariant discrete-time 
systems. This now canonical matrix is the bilinear 
transform of the Schwarz matrix. Liapunov functions 
for linear time-invariant discrete systems can be quite 
simply constructed through tho use of the proposed can- 
onical matrix whose elements, as in the case of the 
Schwarz matrix for the continuous-time systems, contain 
the stability information of the system, attempt is 
made to establish a link betv/oen this canonical matrix 
and tho Jury table test. Several methods of transform- 
ing a given discreto-timo system to one in this canoni- 
cal form are also discussed. 

■Several unsolved problems are suggested for 
further research 



CHAPTER - I 
INTRODtXJTIOH 
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Canonical forms for -system description are used 
extensively in both analysis a,nd synthesis of control 
systems. The significance of these forms is fairly obvious. 

Typically in a synthesis problem, the number of parameters 

% 

to be designed becomes less (minimal) with the use of 

canonical forms. Analysis of systems with such description 

becomes simpler, as elegant analytical results can be 

obtained in such cases. To quote examples, in control 

« 

system literature, the phase-variable and the Jordan 
canonical forms are well known. 

The Jordan form displays the system eigenvalues 
explicitly and hence is ideally suited for the analysis 
problem. However, given a system in any arbitrary form, 
transforming this to the Jordan form, involves the computation 
of the system elgGnva3.ues and eigenvectors and thus becomes 
computationally difficult especially when the system size 
is large. Eiorthor, the problem of degeneracy with repeated 
eigenvalues makes such a trarsf or motion still more difficult. 
The phase -.variable form, on the other hand, is quite simple 
and the corresp-r'^l nc transformetion is comparatively easy. 
But, the elements of the phase-variable form do not have 
as much information (for instance, stability of the system) 

8«! those of the Jordan form. This investigation cohoerhs 
with £ ofinon^<'al furm which has both simplicity In terras 


k 
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of the associated transformation and "enough” information 
about the system behaviour "stored” in its elements. 

Schwarz^l"* in 1956 introduced, for stability 
investigation of linear time-invariant continuous-time 
systems, a new canonical form which was later popularized 
by Kalman and Bertram [s1 and which is now commonly known 
as the Schwarz canonical form. This canonical form has 
been used to construct Liapunov functions [2^ , to prove 
the Routh-Hurwltz stability criterion through the second 
method of Liapunov , to derive optimum transfer 

functions [[4] , to evaluate system performance measures ^,5 ,6 
and to solve the Liapunov matrix equation [j0,9,lO,llJ . 

The Schwarz canonical form for linear systems is 
of central interest in this investigation. Three different- 
but simple r- forms of the Schwarz canonical matrix are 
identified and the similarity transformations -between 
these forms are established. G-iven a system i.n phase- 
variable form, transforming this to the Schwarz form has 
been achieved quite recently j^l3,14j , This transformation 
process has been further simplified in this investigation. 
When the given system is in any arbitrary form, the 
associates brans f* or mat ion to the Schwarz form is described, 

A Tandermondo~like transformation matrix which takes the 
system description from Schwarz to Jordan form has been 
found out, Schwarz canonical form for multivariable systems 
has becii proposed. This resiilt has been obtained bv 
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extending the recent work of Luenberger {jLS j i A simple 
and interesting network interpretation is given to the 
Schwarz canonical form both for single-variable and 
multivariable systems. It is hoped that this hett/^ohk 
interpretation would facilitate both analysis and synthesis 
problems. Added to these theoretical investigations, an 
important practical use for the Schwarz form in 3implli‘ying 
large dynamic systems has been suggested and the simplification 
procedure has been investigated in detail* Several examples 
have been considered for illustrating the method. A new 
canonical form for linear time-invariant discrete-time 
systems has been proposed. This canonical form has been 
obtained thro\:igh the bilinear transformation of the Schwarz 
form. It has been shown to be useful in constructing 
Liapunov functions. The properties of thih oahonical form 
have been thoroughly investigated. A number of unsolved 
problems are suggested forfxnrther research. 



CHAPTER - II 

scmmz CANONICAL form 
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2 a lETRODCOTION 

In this chapter, Schwarz canonical form for 
linear time-invariant continuoUs-time systfenia la Intl^oduced. 
Some of the results regarding this canonical form ■WhioJi are 
already available in the literature are stated* 

In section 2, the Schwarz canonical matrix is 
developed for systems having characteristic polynomials 
with real coefficients. Three different (but similar) forms 
of the Schwarz matrix and the transformations between them 
are also given in this section. Section 3 discusses 
various properties and uses of the Schwarz canonical 
matrix. 

2.B DEVELOPMENT OF TEE SCHWARZ CANONICAL MATRIX [J 
Consider a system described by 

X = AX (2.2,1) 

Let the characteristic equation of (2.2.1) be given by: 

P{A.) = + ii*. + ^n “ ^ (2.2.2) 

The polynomial 

Q(A) = a, + .... , (2,2.3) 

1 3 

is known as the alterant of P( A) in (2.2<l) [l6j[ . 

The test fraction of P(A ) is defined as: 
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P( /\) 


+ 1 + 


^2- "• 


^ V + — 

2 i* 4A + % 


(2.2.4) 


Where 


r 

m 


!ff . -i . } ,.., , .. ( i a il.L > (ni»l,2,...,n) (2.2,5) 


j odd 


1 . 2 ©■'Ten with a » 1 

J ^ 


j even 


( 2 . 2 . 6 ) 


Ij j odd J 

b • t) . ( " 0,1,2, .,,,n ) o ) 

Suhstitutii4 in ^ U-a:) C-^ ^ 


= "bn » 


^k-1 


°n-k+i » 

(k = 2,3,.,.,n) (2.2,7) 


results in 

^l\ ) 

P( A ) 


-b -A + 

n 


(2.2.8) 


4 



to I 
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By applying a theorem in continued fraction 160 


equation (8.2.8) yields 


P( A ) = 


^ n-1 


. , , 0 0 i 


^ . 0 0 


... 0 0 


. (S.2.9) 


i » 0 ^ ji 


where 


== |§^ - A I 


\-i ° 


... “xA h 2 ^ 

... —1 - ^ I 


0 , 0 

0 . .. 0 


^n-3 • • * ^ 


(2.2.10) 


(2.2ai) 


... 0 
... -1 


2 in (2.2.11) is similar to A in (2.2.1) and it is 

one of the three forms of the Schwarz matrix. The other 

two forms are denoted "htr q onri « t i. 

X« ueno-cea by and Sg and they are given 


as : 
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-b --b 
jci-l n 


(2. 2. IS) 


-2 


(2.2.13) 


— J 

Observe that Sg and Sg are given In | ij whereas 

Kalman and Bertram UJ ha™ made use of the form s,. 

In this report also, §j is mostly used and the litter 

B is interchangeably used to denote s^. Whenever we 

refer to the Schwarz matrix without any qualification, we 
inean Sj, only. 

The similarity transformations between S and 
S 3 and between S^ and will now be oonsilered! Let 

§2 = T;' 


(2.8.14) 
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S-1 

-5 - 1-3 


(2,2*15) 


Tiien 

‘^3>is 

(To) 

21 

‘Es’ii 

'SsV 

f5s'i(i+i) 


‘-s’(i+i)i ' '‘-s’ki+D ^ 


Tg as given in (2.2.16) is tridiagonal. Except for 
the element (To) , To is skew- symmetric. For n=4, 

JlIiI — <0 


0, i=l ,2, ... ,n-l 
n 

- J 'x. (-h, ) 

k=2 ^ k^ 




f i=2 ,3, . . . ,n-l 


(2.2.16) 


Se = 


(-b^) 


0 (-bg) 




0 0 (-bgX-bg) -(- bgX-bg) (-b^) 


The transformation matrix is a cross-diagonal matrix 

(i.e., the only nonzero elements lie on the cross- 
diagonal) and is given by: ^ 

(Js),, ' = 0, ,3 / n-1+1 i 


(S3)lj ' 

‘^s’Kn-l+l) 


= (-1)“"^*^, (1=1,2,. ,.,n) 


(2.2.1V) 
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For n = 4, 


?3 = 


0 

0 

0 

-1 


0 0 1 

0-10 
10 0 
0 0 0 


In cJiapter IV, the concept of network interpretation to 
the Schwarz forjtn is effectively -used to derive a simpler 

between and Sg. 

3.3 PROPERTIES AW IBES OF THE SCHWARZ CAHONICAl FORM 
Consider a systeia described by 



y “ B y + f u 

where B is as in (3.3,12) and 


(3.3.1) 


f = ( 0,0,.. .,0,1) 


T 


(3.3.2) 


Elaborately writing the differential equations in (2,3,l) 


y- 


= y? 


-bjyi + Yg 

■'’sS's * 

->> 3^3 + 


n-3 

-\-4 

^n-4 

4 

^-2 

' = 

-■b a 

y 

4 

y 

n-2 

n-3 

n-3 


■^n-1 

n-1 = 


yn-P 

4 

yn 




(2.3.3) 
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Tiiese diff erontial equations can be equivalently represented 
by ftieans of a simulation diagram as in Figure (2,3.1) (In 
bills Figure, n is taken as an even number, A similar 
diagram may be dravm for n odd case also). Except for 
tile first and tile last differential equations in (2,3,3), 
otJiers have the same pattern. The simulation diagram 
also exhioits such a pattern between the dotted lines i_See 
Figure (2*3.1) i.e., between the first and last integrators. 

Such simulation diagram may be useful for computing the 
response of a system in Schwarz form by using analog computer 
or digital computer (For instance, "Pactolus” program may 
bo used). In fact, this has been used for the example 
systems in chapter T, 

By modifying the results in cohtihUed fraction 
method f_ 16 l , Schwarz [l] has obtained the following: 

The system in (2.2,1) is asymptotically stable if and only 
if all tho t^’s (i==l,2,, ., ,n) in (2,2.11) or (2.2,12) or 
( 2.2.13) are posits ye^ The number of eigenvalues of A 
in (2,2,1) with negative real parts is given by the number 
of positive terms in the sequence b^, bj^bj^_^, ’ * * * 

"'^n'^n-l'^n- 2* * *'*^1 * is an important result and for 

future reference, let us call this as Schwarz theorem. 

Kalman and Bertram {_2] have proved this theorem through 
the second method of Liapunov. Consider a system 

y = B y 


(2.3.4) 














'■nl ^ts Liapunov function 

-rp 

Y = y L y 

v/here 


IS 


(3.3.5) 


(2,3.6) 


Then 


V 


2 2 
-2b „ y 
n ''n 


(2*3.7) 


which is negative semidefinite . Also Y cannot be 
Identioally zero unlees = 0 (1=1 ,2, . . . ,n) [ 4 ] . Thus 
from the second method of Liapunov, the system in (2.3,4) 
is asymptotically stable if and only if L in (2.3,6) is 
positive definite. Hence Schwarz theorem is proved. 

Parks IJ5,4:J has shown that the Schwarz elements 
bp's (i=l ,2 , . . . ,n) are related to the Hurwitz determinants 
/h^’s (i=l ,2, , . . ,n) as under: 


h = 


/> 


n-i + 1 


A 


n-i-2 


, (i=l,2, . . . ,n) 


(2.3.8) 


/S A n-i-l 

with k£.l. It may be easily verified that 

r in (2.2.6) is expressed in terms of the elements 
m 

in tho first column of the Routh array, R . ( j=l ,2 , , , . ,n«*-l) , 

Ij^ 

of (2,2.2) as follows: 


R 


= 

m 


ml 


R 


(m+l)l 


m » 1,2, ... ,n 


(2.3.9) 
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with = 1. From (2.2.7) and (2.3.9) it is not difficult 

to see that 

V ^( n--iH-2)l ^ 

i jj ’ i — l;2j...,n ^2.3.10) 

(n-i)l 

with = 1 when k ^ 1. From (2*3.8), (2.3.10) and 

Schwarz theorem or the result of the second method of 
Liapunov (i.e., Ii> £) > it is inferred that all the first 
column elements of the Routh array or all the Hurwitz 
determinants must he positive. This is how Parks L5,4j has 
established the link between the second method of Liapunov 
and the Routh-Hurwitz stability criterion. He has also 
used the Schwarz canonical form to derive ’optimum transfer 
functions’. For details, refer to jjrJ. 

From equations (2.3.5), (2.3,6) and (2.3,7) it is 
clear that the Schwarz canonical form can be used for 
constructing Liapunov function for a •' system* Suppose 
the given system has description as in (2.2,1), Apply a 
transformation 

y = P X (2.3.11) 


* A in (2.2,1) may be in any arbitrary form or in 
a special case in phase-variable form, According to the 
notation used in latter chapters, P is the transformation 
matrix such that B = P A P“1 where A is in phase -variable 
form. Instead if a Ts^in any arbitrary fofm, H is the 
transf ormation matrTx such that B = ^ H, Hnwever, in 

(2,3.11) we just use P without assuming^tEe form of A, 



Tiion tile Liapunov function for tiie system in (S.2.1) can 
be easily written as 


1 t 1 (2.3.12) 

in view of {2.3,5} and (2.3,11). Tlie transformation 
(2.3.11) will be discussed elaborately in chapter III, 

Tiio Schwarz canonical form can be used to evaluate 
the system performance measures of the type to be given 

■7 

below L5>6 . But this requires that the system description 

is to be in Schwarz form [See equation (2.3.4)J , Assume 

that the system is asymptotically stable. Consider now the 

* 

performance measure 

OQ 2 

I = 7 y dt (2.3.13) 

o ^ 

Integrating (2.5.7) between o and co With respect to 
time results in 


I = 




( Y(od ) - T{o) ) 


(2,3.14) 


Noting that for an asymptotically stable system, y(oo) = 0 
and substituting (2.3,5) in (2,3.14) 

I = y! y'^(o)Ly(o) (2.3.15) 


n 


2b‘ 


n 


In roforence fsj , this method has been extended to 
evaluate performance measxires of the type 
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Recently . 7 ,8 ,9 ,10 ,11 ,1S_, there has "been a great 
deal of interest shown towards solving the Liapunov 
matrix equation 

K + K ^ = -G (2.3.17) 

whore K>£ and G^O. The solution of the equation (2.3.17) 
may be useful for calculating functionals of motion [l7j , 
estimating overshoots and settling times [isj and deriving 
pseudo-optimal control policies for tho vector u which 
will roturn the system 

i 2 (2.3.18) 

to equilibrium as quickly as possible following an Initial 
distiirbando 1 18] * Schwarz canonical form has been used 
to determine the solution of (2,3.17) quite simply. For 
details, refer to [8,9,10,ll7 . 

2.4 COROLIBION 

Some of tho uses for the Schwarz canonical form 
disGiissed in this chapter rely upon our ability to compute 
the transformation matrix P [^Seo equation (2,3,11)J . This 
transformation is discussed elaborately in tho next chapter. 
In chapter 7, it is shown that the Schwarz canonical form 
can be used to simplify largo dynamic systems. A possible 
topic for fiirthcr investigation is the modeling of the 
physical systems to get the system description directly 
in Schwarz form. 
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CHiiPTM - III 

TRAmFORMATIONS ASSOCIaTID V/ITH SCHWARZ CANONICAL FORM 
3.1 INTRODUCTION 

In the lost chapter, some of the uses for the 
Schwarz canonical form were discussed. However, it was 
not stated as to how a system in any arbitrary form could 
be transformed to the Schwarz canonical form. This problem 
is treated in this chapter. Schwarz |_ll used a series of 
Giementary transformations to transform any given matrix 
to its Schwarz form. However, no simple similarity 
transformation like the one between the arbitrary form 
and the phase-variable form [l9j was reported till recently. 
Recently there has been a great deal of interest shown 
towards finding such a transformation to the Schwarz 
canonical form [l3)14,S0|Sl3 . 

The contents of this chapter will how be briefly 
stated with the help of -fche Figure (3»lil)4 Section 2 
discusses the transformation between the phase-variable 
form C and tho Schwarz form B i.e., the link No. 5 
in Figure (3.1,1) . The elements of tho corresponding 

transformation matrices P and can be explicitly 

expressed in terms of either the oloments of the 

Routh array or the elements b^»s of B. A critical 
discussion on these methods of determining P and P"^ 
is also contained in section 2. Section 3 deals with 
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tiie transformation from the arbitrary form A to the 
Schwcirz form B (link No.4K The transformation between 
the Schwarz and the Jordan forms is the subject of 
section 4 (link No. 3). An elegant method of determining 
N (link No.l) is given in L2o} , Of course, the 
transformation in link No. 2 is achieved through the 
familiar Vandermonde matrix, V. In section 5, Schwarz 
canonical form for multivariable systems is proposed and 
the transformation to this canonical form is also discussed. 

3.2 THE TRAI^FORMATION EROM THE PH/i^E- VARIABLE FORM 
TO THE SCHWARZ FORM 

Consider a system described by 

X = C X (5.2.1) 

where C_ is in phase-variable form. It is desired to 

find a transformation 

y = P X (3.2.E) 

which takes the system description from (3,£.l) to 

J = B y (3.2.3) 

whore B is in Schwarz form. In other words, 

B “ F C (3,2.4) 

I*b has heeh found ri3,l4j that the elements of P and 
P"^ can be expressed Ih tehms of the elements of the 
Routh array, R. .’s, of the system. The resiilts will 
now be stated. The matrices P and P”^ have the 
following interesting properties: 
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1. They ore lower triangiilar matrices. That is, ail the 
oloments above the main diagonal of P and are 

zeros . 

3. All the elements on the main diagonal of P and P 
are equal to 1, 

3. All the odd-ordered, subdiagonals of P and P“^ 
contain zeros. 

4. The even-ordered subdiagonals of P and P”^ contain 

nonzero elements and these elements are determined 
from the elements of the Routh array of the 

X J 

system, 

5. If P^ is the m * n transformation matrix P, 

then , k = l,2,...,n-l (That is, Pj^ is the 

transformation matrix for a system of order k) are 

the submatrices located at the lower right corner 

of P . Such submatrices will be hereafter referred 
-n 

to as tho ascending principal minors of P^ of 
appropriate order. The matrix P~^ also has this 
property. 

Determination of the elements on the even-ordered 
subdiagonals of P and will be discussed next. To 

begin with, the results obtained in references [^isj and 
[14] will be stated briefly. Then |t will be shown 
how the samo results can be put in a simpler and more 
useful form. 

Consider first the matrix P. For convenience, 
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let n = 8. The P raatrix for this case is shown below: 


P = 


1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

^72 

R 

71 

0 

1 

0 

0 

0 

0 

0 

0 

^62 

0 

1 

0 

0 

0 

0 


R 

61 







^53 

51 

0 

^58 

R5I 

0 

1 

0 

0 

0 

0 

^43 

R 

41 

0 

R.r, 

42 

R 

41 

0 

1 

0 

0 

^34 


^33 


^32 




^31 

0 

^31 

0 

R 

31 

0 

1 

0 

0 


0 

^23 

0 

!!! 

0 

1 




^21 


^21 


1 


(3.2.5) 


Observe from (3.2,5) that the elements of P on the 

second subdiagonal arc computed ]3y dividing the corresponding 

oloraonts in the second column by the elements in the first 

column of the Routh array, in fact, the last element in 

^22 

this second subdiagonal is , the element just before this 
is etc. and the fi:^^ eleftient in this subdiagonal 

^31 



( in this case it is 
computing the fourth 


IS 


^(n-l)8 

^{n-l)l 


R 

R 


72 

71 


SI 

) , Similarly, for 


suhdiagonal elements 
of P , Routh array elements in the third and first columns 
are to he considered and for the k-th suhdiagonal ( k even ) 
elements of P , Routh array columns ( k/S + 1 ) and 1 
are to he considered. 

Consider nezt P“^, The elements on the second 
subdiagonal of P“^ are ( from the lower right corner to 
the upper left corner ) given as. 


R 


as 



R 

42 


R 

{n-l)2 


^21 ^31 \l ^(n-l)l 

Note that these are same as the elements on the second 
subdiagonal of P with a negative sign. From the 
lower right corner to the upper left corner, the elements 
on the foincfch suhdiagonal of P~^ are the 2x2 
determinants 


^42 


^58 


^62 


^(n-l)2 

. — ... 1 


1 


1 


1 

^41 


^51 


^61 


^{n-Dl 


5 

R R 

5 

R R 

J • • • > 

R , \ R / 

23 22 

1 

33 32 


43 42 


(n-3)3 (n-3)2 

R__ R_ . 


R R 


R R 


R / K R / \ 

21 21 


31 31 

' ! 

41 41 


""(n-S)! ^(n-3)l: 
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■which are illustrated in the following; 



The elejnents on the sixth subdiagonal of P“^ are 3x3 


determinants which can be obtained from P in a similar, 
manner L143 , Note that this procedure of determining P 
involves com.putation of determinants of order upto[^^2 
where the square brackets refer to the integer operation, 
i.e,, [Jx3 means the largest integer contained in x. 

These resiiLts will now be put in a more compact and 




23 


Note that in (3.2,6) P is represented in a quite simple 
and compact form and that p/ ^ , , s = l,2,...,n-i in 
(3.2,6) are the elements along the i-th subdiagonal of P. 
Consider next P"^. For convenience, let 


Then 


q * 

ij i n X n 


j > i 


i odd 


[3,2.7) 


■(s+i) s 


= ^ 




+2-m 


m=l (n-s-i+2)l 


^(2m+s-2)s 


i even 


i = 1,2, . ,n-l 

s=l,2,,,,,n-i ^ 

The following comments are now in order. First, (3.2,7) 
gives a compact representation of ^ . Second, ag before 
'^{s+i)s ’ ^ ~ l,2,,,.,n-i are the elements on the i-th 
subdiagonal of P“^ . Third, computation of determinants 
is avoided in this representation. The computational 
difficulty is now reduced to a large extent because in 
the process of computation, the knowledge of the previously 
computed elements is used effectively [See ecuatlon(3.2,7)_J . 
Equation (3.2.7) is especially useful for machine computation. 
In the rest of this section, it is aimed to determine 
P and p~l in a still simpler way. First, the elements of 
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1 

P and P“-^ will "be expressed in terms of the Schwarz 
elements b^’s (i=l , 2 , . . . ,n) . 

Consider the elements p./. ,i=3,4,...,n of 

the second lower diagonal of P, They are given by 

Pi(i_E) = ^ (i=3,4,...,n) (3.2.8) 

(n-i+ 2 )l 

It is quite simple to see that 
^(n-l )2 " ^(n+l)l 

R u TD (3.2,9) 

(n-i+2)2= ^(n-i+3)2 ^(n-i+2)l^^ 

^(n-i+4) 1 
R 

(n-i+3)l 

( i“4 j 3 } . * • j n) 

Referring to chapter II, 


with Rj^j^ = 
(3.2.10) , 

^i(i- 2 ) 


^(n-i 48 )l 

R 

(n-i ) 1 

1 when k ^ 1 


i ~2 



From (3.2,8), 


i = 3,4, ... ,n 


(3.2.10) 
(3.2,9) and 

(3.2.11) 


Consider now the fourth subdiagonal of P. The elements 
along this diagonal, Pi{j[_ 4 ) > i = 5,6,...,n are given 
by: 

R 

^ .In -ljll g ^ (1 = 5 6,....n) (3.2.13) 

i(l-4) R, . ... 

(n-i+4)l 


Using the relation between the Routh array elements, 



S5 




^(i-l)l 


(3.3.13) 


it is very simple to show, from (3.2,9) and (3,2.10), that 


b p 
3 ^31 


51(1-4) 


5(1-3)(1-4) ^(i-4) (JL.5) 

(i=6 ,7, . . , ,n) 


(3.2,14) 


This prooedTjre can be repeated and in general 
^0 j >i 

0 i+J odd 


^(i-2)l ^(i-2) > I i+j 


(i-2)j ^(i-2) ^(i-l)(j-l) 

• i 3 , i 1 

... y 

{ i»d ~ lj2,,,,,n ) 


even 


(3,2.15) 


In words, any element p. is found by multiplying the 

ij 

element situated in the same column and (1-2) -th row by 
^i 2 adding it to the element situated in (i-l)-th 

row and (j-l)-th column. Of course, the first row is to 
be first filled in as (1, 0, 0, , 0, 0) . 

It may be similarly shown that the elements 
of P“^ are given as: 
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0 j > i 

0 j + i odd 

1 i = j 


"^1 ^(i-l)2 » ^ ^ 


3 = 1 


> (3.2.16) 


'^(i-l)(j-l) ■*' '^{i-l){j+l) 


i + 3 
even 


I ¥ 1 , i i 

In words, first write down the first row of P**^ as 
( 1 ,0,0 , , . . ,0,0) , Then any element is found by 

multiplying the element in the {i-l)-th row and (j+l)~th 
column by (-b.) and adding it to the element in the 

J 

(i-l)-th row and the (3-l)-th column. Power £2oJ has 
suggested a flow graph rule to compute P”^. 

It is quite interesting to observe the following 
points. First, if Pn(PQ^) is known as in (3.2,15) 

[(3.2, 16) J , then Pj^(PJ^^) where k^n is the descending 
principal minor of order k of (P^ ) , 

Second, if P is given as in (3,2.17) where 
are to be expressed in terms of hj^’s as in (3.2,15), then 
the elements of the second subdiagonal of P 


~1 


from 


the upper left corner to the lower right corner, are: 


-P 


31 


'i^42 


“P 


53 


j ♦ * 


’ ■^n(n-2) 


The fourth subdiagonal elements, from the upper left 
corner to the lower right corner, are given by the 2x2 
determinants 



E7 


' ^31 

1 


®4S '■ 

9 

P53 1 

i 


P 1 

(n~2) (n-4) 

I’SI 

1 

^53 

1 

®62 ®64 


^73 ^75; 

i 

^’n{n~4) Pn{nr2) 


as illustrated in (3.E.17), 



... (3.3,17) 


TJiese results are exactly analogous to the previous ones 
except that now everything starts from the upper left 
corner. This is because of the way in which the bj^*s and 
the are related , 

Finally, it is observed that only the first column 
of the Routh array is used to compute the bj^'s from 
(3,2,10) as against the previous method which uses all 
the columns of the Routh array, Eq.ua tions (3,2.15) and 
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(3.2.16) are exceedingly simple to remember and use in 

order to determine P and respectively, ismmber of 

multiplications involved in the computation is considerably 
reduced. It is for these reasons, equations (3,2,15) and 

(3.2.16) are recommended for computing P and 
respectively. 

one more method to compute the transformation 
matrices is available [20J , After computing the elements 
bj^’s of B ( from the Routh array ), apply Power’s 
transformation L2o3 to B so as to transform B to its 


phase- 

-variable form. 

Thus 



0 = 

¥ 

B 

(3.2.18) 

where 

M can be 

found by Power’s method. 

Therefore , 


p-1 


M 


and 

P 

ss 


(3.2.19) 


As this method is not superior to the one discussed above 
^equations (3,2,15) and (3,2,16) J, detailed presentation 
of this method is not aimed at. However, further details 
may be obtained from [boJ . 

3.3 IRANSFORI^TION FROM ARBITRARY FORM TO SCHWARZ FORM 
Recently there have been many papers [19J discussing 
the transformation of a controllable system whose description 
is in arbitrary form to the phase-variable form. This 
section, following Tuel’s [B2j approach* derives A ' 
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transformation wMcii takes a controllable system in 
arbitrary form to one in Sciiwarz canonical form [21 J , 

Tile transformation is extended to unforced systems, 
thereby pointing out the difference between forced and 
unforced systems from the viewpoint of such a transformation. 
A particular case wherein the given system has the phase-' - 
variable form description is considered and under a certain 
condition, the corresponding transformation is shown to be 
the one discussed in the last section. 

Consider a controllable system 

X = A X + u (3.3. 1)_^ 

where A is some arbitrary n x n matrix. A transformation 

X « H y (5.3.2) 

is sought such that (3,3.1) becomes 

y » B y + f u (3.3.3) 


where B is ih Schwarz form and 

f = ( 0, 0, . . . , 0, 1 )'^ 

From (3.3,1), (3,3.8) and (3,3,3) it is seen that 

H B = A H 


(3.3.4) 

(3.3.5) 


and 


let 


H f = t 


§ = r 1 



(3.3.6) 

(3.3.7) 


Then (3,3,5) yields 
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\ 

o 

o 

ft ft 

. 0 0 



-b^ 0 1 

ft *> 

. 0 0 



0 -bg 0 

ft ft 

. 0 0 


fei; -2' “* 

— , » 

« # ft 

ft ft 

••ft « 



ft ft ft 

ft * 

ft ft ft 



0 0 0 

ft ft 

0 1 



0 0 0 
L_ 

ft ft 

“^n-l"^n 




ft ft ft 

(3*,3*8) 

i .e . , 





k bi^A h ! .. . 'u li 7 = 
L J-,- -S, -nj 

= [^-b^ 

r 

^ h ^ 

^2 i-'S >• • • 

;-n-2 “^n-1 ^n! 

i 


-n-l -^n ^?-nl 

ft ft ft 

(3.3.9) 

From (3.3.6) it is clear that 




-n 

= b 



(3.3.10) 

Equations (3.3.9) and 

(3,3,10) give 

the following results: 

h 

= b 




-n 





h„ T 

-n-l 

~ - ^ ^n 

b 


(3.3.11) 

h , 
-n-k 

= il h , , 

- -n-k+1 

+ b 

n-k+l 

-n-k+E 



(k— 2 > 3 , , , , 

,h-l) 




Tile controllability condition imposod on tbe pair (ii , b) 
implies that tlie composite matrix 

- = f- ' ^ 'b 1 • . . hi (3.3.1; 

^ < / " 1 

is nonsingxilar. From (3.3.7) and (3.3.11), H is written 
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H = [4^“^ t+lD^ b+ ... ! ... -'- -‘*’^n-'-l 

I } i ' 

... (3.3.13) 

By performing elementary column operations on H , the 
equivalence between H and Q can be established. Hence 
F is nonsingular if and only if the aforementioned 
controllability condition is satisfied. 

It is interesting to note that the knowledge of 
the bj^'s j i = l,2,..,,n, the elements of the Schwarz 
matrix B, is required to determine H in (3*3*13), Tuel's 
transformation matrix |_2^ K ( K transforms the system 
(3.3.1) to the one given by (3*3.3) where now B is in 
phase-variable form ) also involves the elements of the 
characteristic equation of a* a.n alternative procedure 
to construct K without aprlori knowledge of the 
coefficients of the characteristic equation of a is 
given by the Luenbergor [isj , Such a procedure does not 
seem possible in the case of H matrix because of the 
fact that the bj_’s in B matrix are distributed 
throughout its rows and columns. 

From the main assumption, the system (3,3,3) is 
controllable. This fact leads to the following conclusion: 
**. sufficient condition for a system which is in Schwarz 
form to be controllable is that the system input vector 
must have a 1 in the n-th position and zeros everywhere 


else. 
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Consider now a free system 

X = A X (3.3.14) 

Seek a transf ormation 

X = H y (3.3.15) 

such that 

1=^1 (3.3.16) 

H in (3.3.15) can now he determined as follows: Choose a 
vector h such that the pair (A , h) is controllahle. 

Then form H as given in (3.3.7) and (5.3,11) with this b. 

Considering the forced system, once the forcing 
vector b is specified, the transformation discussed in 
this section is possible if and only if the pair (A , b) is 
controllable. But for free systems b can be so chosen 
as to ensure the controllability of the pair (^^ > b ) . In 
fact, there may be many transformation matrices H 
transforming (3,3,14) to (3,3,16), 

Suppose A in the system (3,3.14) is in phase- 
variable form (denoted by G) . In section 2, a transformation 

y = P X (3.3.17) 

such that 

P C P~^ = B (3.3.18) 

was discussed. The elements of P in (3,3,17) are 
determined from the knowledge of the elements of the South 
array of (3.3.14). Whereas the choice of P is unique, 
that of H is not so owing to the freedom available in 
choosing b. It is intuitively felt that a particular choice 
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of b may make P"^ to be identically equal to H, 

It can be establislied that when b == f , P“^ = H. It 
is quite simple to show that the system (3.3.1) with 
^ = 0 and b = f is controllable. Determination of H 
with " = £ and b = f can be accomplished from the 
knowledge of either the Routh array elements or the b^'s 
(See section 3,2). This can be used to determine 
in terms of the elements of the Routh array. In fact, 
this is how the author obtained result in (3.2.7). 

A 

■fi.n example will now be considered. 

Example: 

Consider 


• 






"^1 


0 10 0 


^1 

^2 


0 0 10 


^2 


xs 




^3 


0 

0 

0 


^3 

X4 


-3 _1E -6 -4 

j 

1 1 

^4 





The characteristic equation for this is 

+ 6>^4l2>+3 = 0 

Routh array is then formed, 

16 3 

4 12 

3 3 

8 
3 


TJsing the relationship between the Routh array elements 
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(R. .’s) and th.e Schwarz elements 
1 J 


P = 




'‘si 


= 

4 





^3 

“31 


= 

3 





^2 = 

R /R 

41^ 21 



2 





"l = 

“3i/“31 


= 

1 




15) 









1 1 

0 

0 0 


1 


0 

0 

0 

0 

1 

0 0 


0 


1 

0 

0 


0 

1 0 


1 


0 

1 

0 

0 

(b^+bg) 

0 1 


0 


3 

0 

1 


From (3.S.16) 


1 

0 

0 

0 


1 

0 

0 

0 

0 

1 

0 

0 


0 

1 

0 

0 

-bl 

0 

1 

0 


-1 

0 

1 

0 

0 

-(b^+bg) 

0 

1 


0 

-3 

0 

1 


To compute H , let h = f. Thus, using (3,3,11) 


-4 

f 

= ( 

0 

0 

0 1 ) 

-3 

A 

J + f . ( 

0 

0 

1 0 ) 


A 

&3 + *> 354 = ( 

0 

1 

0 - 3 ) 

^1 = 

A 

62 ♦ Vs" < 

1 

0 

-1 0) 

Hence 







10 0 0 

0 10 0 

-10 10 
0-301 
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Observe that H = 

3,4 TRANSFOro^TION FROM SCHWiiRZ TO JORDAN FORM \z^'] 

This section pres(§ints a transforiaation which takes 
a system from the Schwahz Tohm to the Jordan form when the 
system has (i) distinct eigenvalues, (ii) repeated 
eigenvalues having simple degeneracy. 

Consider a system described by 

2 = B X (3,4.1) 

It is desired to find a transformation 

_z = T z (3.4,2) 

such that 

z = J z 

where J = T“^ B T is the Jordan canonical form 
corresponding to the Schwarz form B. Referring to 
Figure (3.1.1) it can be readily seen that 

T == P V (3.4,3) 

In section 3.2, it was stated that the elements of P 
can be easily determined in terms of b^’s by a simple 
rule fsee equation {3.8.15)J * However, it is also 
possible to express the elements of P in terms of the 
b^^’s explicitly [^Note that relations in (3.2.15) are 
algorithmic in nature J , 
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ij 


p = 

0 
0 
1 
0 




ij 


J > i 


i + j odd 
i = 3 

j + 2 i _ 2 


r (3.4.4) 


( . • » ) 


k=l f=k+2 

i ;} 

(i» j = 1, 2, 3, n) 




i*3 

even! 


This simplification helps in expressing T solely in terms 
of the elements of B and its eigenvalues. 

Distinct Eigenvalues 
If T is represented as 


~ " [■‘^ij ] > {3.4.5} 

then by performing the matrix multiplication T = p v , 

it can be shown, in the case of distinct eigenvalues, 
that 


t — \ I~1 . 

^13 ■ t 


i-2 

k=l k 


vi-3 

^3 * 


i-4 i~2 


7 


b, b, 






1-5 




i-6 



i-4 i-2 


s z: 

-t=k+2 m=f+2 





( IjD ~ 1,2, ...^n) 


(3,4,6) 


where \ = 0 if k<0 and denotes the i-th 
eigenvalue of the system in (3.4.1). An algorithm for 



generating the is given below: 

X J 
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hk 

= 1, 

(k = 1 ,8 , . . . ,n) 

"^Bk 

= Aj, 

, (k = 1 ,2 , . . , ,n) 

t 

Ik 

= /'k 

^(i-l)k ^i-2 ^(i 


(3.4.7) 


(i = 3,4, ... ,n) 

(k = 1 ,2 , . . . ,n.) 

Typically in the case of a fifth order system, if 
T = ^t ^t [t ^t|t 

and ^ (k = 1,2,.,. ,5) are the eigenvalues, then 


1 

Ak 

^ 2 

Ak ”1 

* <'=1 ‘'2 * ’> 3 >>k + Vs 

( k = 1,2, ... ,5 ) 


Case II : Repeated Eigenvalues with Simple Degeneracy 
lor simplicity, let us assume that B has one 
repeated eigenvalue with multiplicity m, other 

eigenvalues being distinct. Let 


T 



t 

-m 


-m**"! I 


3 


In forming the matrix product T = P V » J. is the 
generalized Tandermonde matrix [^2oJ . The resulting 
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transformation matrix T is such, that the columns t^ , 
-m+1 * -m+2 ’ ’ ~n computed using (3.4,7) and 

1 d^~^ 

^ t^^ ^ 2,3, ...,m) (3.4,9) 

(i-l)I dA“ 

The same procedure is valid even when B has more than 
ono repeated oigonvalue, provided each oigenvaluo has 
simple degeneracy. 

Example : Consider the system (3,4.1) with 

0 10 0 
10 10 

B = 

0 0 0 1 

0 0-2-3 

Here b, = -1 , h. = 0 , = 2 and b^ = 3. B has 

eigenvalues -1, -1, 1 and -2. It can be checked that 
the eigenvalue -1 which has a multiplicity m = 2 has 
a simple degeneracy. By using the procedure outlined 
above, T can be computed as: 



Now it can be easily verified that T“^ §. 2 “ £ where 

-110 0 
0-100 

J = 

0 0 10 

0 0 0 -2 
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It may be noted that the transformation T given 
above has some resemblance to the Vandermonde matrix. 

3.5 sommz CANomcAL mm foh mtivariable sYsm© 

In this section, Schwarz canonical form for 
multivariable systems is suggested. A procedure for 
transforming the system description from arbitrary form 
to the Schwarz form is outlined. An example system is 
considered for illustration. 

Luenberger [15 j has recently proposed a transformation 
which takes a controllable miiltlvariable system in arbitrary 
form to a set of coupled subsystems each of which being in 
phase-variable form. In this section, the transformation 
( P ) discussed in section 3,8 is applied to each such 
subsystem to bring it to the Schwarz form. Combining 
these two transformation processes , a single transformation 
may be obtained; this, therefore, takes a controllable 
multivariable system to a set of coupled subsystems each 
of which being in Schwarz form. 

Let the multivariable system be given by 

X = A X + F u (3.5.1) 

whore a is in arbitrary form, i^pply Luenberger’ s 
transformation [15 j 

X = S X (3.5.8) 

and the transformed system becomes 

^ = A 'x + P u (3,5.3) 
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Tile control input u can also be transformed as 

V = G u (3,5,4) 

where G is an upper triangular matrix so that 

I u = F V (3.5.5) 

where 

0 
0 

• 

• 0 

• 

1 

0 
0 

A 

F = 2*0 

« 

0 
1 


0 

0 

0 

M « 

0 


Suppose (3.5.3) has m coupled subsystems 
corresponding to the submatrices , Ag , 
of jx. Each of these submatrices 

is in phase-variable form. Let , P^ , P^ be 

^ — 1 

the transformation matrices such that P^ P^ = 
k = l,2,..,,m where B. is in Schwarz form. Knowing 
can be computed as discussed in section 3.2. 
Now, form a composite transformation matrix P as 
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?2 


0 


0 



(3.5.7) 


Clearly 

y=PSx=Px (say) (3.5.8) 

transforms (3.5,1) to 

y - B y + P V (3.5.9) 

A A 

Note that in (3.5.9) , the fact that P P = F has been 
used. The multivariable system in (3.5,9) is said to be 
in Schwarz canonical form. In (3.5.9), B has m 
submatrices Bj^ , k = l,2,...,m and each of these 
submatrices is in Schwarz form. Observe also that (3,5.9) 
still has coupling terms between subsystems though they 
might have been modified because of the transformation 

(3.5.8) between (3.5.3) and (3.5,9). 

It is interesting to ask as to whether P in 

(3.5.8) can bo directly determined (i.e, without having 

A./ 

to compute a), it seems to the author that it is not 
possible for (i) the dimensions of the submatrices 
or B are not known apriori (i.e. before computing 

iC 

/\ 

^) and (li) p does require the knowledge of the 
Routh array of the submatrices. 
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Example: Consider the system, in (5.5,1) with 



Luenberger’s transformation matrix S can be found as: 



0 0 0 8 -125 108 J 

The matrix triple product S A results in 

0 1 0 0 0 0 

0 0 1 0 0 0 

-6 -11 -6 0 0 0 

A = 

0 0 0 0 1 0 

0 0 0 0 0 1 


0 


0 


0 -120 -74 -15 
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Also 

0 0 
0 0 

^ 10 
F = S F = 

0 0 
0 0 

L ° ^ 

II 

Note that in this simple example, A does not have any 

A 

coupling terms and F = ? (i.e, G = I). Clearly 



and 


0 10 
= 0 0 1 
-120 -74 -15 

Now and A2 can be transformed to their Schwarz 

forms and Bg respectively and the corresponding 

transformation matrices and Pg can be determined. 

The composite transformation matrix P is obtained as 





1 

0 

0 1 
i 

0 

0 

0 




0 

1 

1 

0 1 

1 

0 

0 

0 


0 


1 

0 

1 

1 1 

0 

0 

0 

-1 I 





1 _ 




— — f 
0 : 

p 


0 

0 

0 I 

1 

0 

0 

- “ ' 

-2j 


0 

0 

1 

0 ' 

0 

1 

0 


0 0 0 


8 0 1 
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The transformation matrix 

A 

P = 

P S 

is 

computed 

next. 



1/6 

-1/2 


1/2 


0 


0 

0 



-1/6 

1 


-3/2 


0 


0 

0 

A. 


0 

-5/2 


5 


0 


0 

0 

P = 


0 

0 


0 


1/2 


-5 

3 



0 

0 


0 


-8 


25 

-18 



0 

0 


0 


12 


-165 

138 

Thus 

y 

A 

= P X 

transforms 

the 

system (3 

.5, 

1) to 

(5.5.9 

where 

A A' 

P = F 

(computed 

above) and 







0 

1 

0 

0 

0 


0 





-1 

0 

1 

0 

0 


0 





0 

-10 

-6 

0 

0 


0 



B = 













0 

0 

0 

0 

1 


0 





0 

0 

0 

-8 

0 


1 





0 

0 

0 

0 

-66 

-15 















3.6 GONCLDSION 


Various transf ormations related to tiie Schwarz 
canonical form have been quite extensively discussed in 
this chapter. Also in the last chapter, transformations 
between the three Schwarz forms were dealt with. However, 
a gap still remains to be filled in and that is, finding 
the inverse of the transformation matrix ^ (between B and • 
Analogous to the result available for the inverse of the 
Vandermonde matrix [ 24 J , it may be possible to device a 
method of finding ‘^so some practical use may be 

found for Schwarz form for multivariable systems. 
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chapter - 17 

MBT'TOK IWTlRPRET*iTION FOR SCHWiiRZ CiiRONICAL FORM 

4.1 INTRODUCTION 

In chapter II, the Schwarz canonical matrix was 
developed through the use of continued fractions. Also, 
the Schwarz matrix is tridiagonal. So it is reasonable 
to expect a ladder network interpretation to ihb Schwarz 
canonical form. In this chapter, the three forms of the 
Schwarz matrix S^ , Sg and Sg introduced in chapter II 
are associated with resistively terminated L-C ladder 
networks. Section 2 derives elaborately such a network 
interpretation. As a byproduct, simple transformations 
between and Sg and between S^ and §3 are 

obtained. The possibilities of using such an interpretation 
for both analysis and synthesis of systems are also 
Investigated in the same section. Section 3 extends 
this network interpretation to multivariable systems in 
Schwarz canonical form. An example is also contained in 
this section. 

4.2 NECTORK INTERPRETATION TO SI NGLE- VARIABLE SYSTMIS 

IN SCHWARZ FORM 

In this section, it is shown that the Schwarz 
matrices , Sg and §3 are associated with resistively 
terminated L-C ladder networks. As a byproduct of such 
a network interpretation, the transformations between > 
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§2 and S3 are derived in a simple way. 

Suppose title transfer function of a system is of 
the form: 


G(s) = (4.S.I) 

a s^ + a, 3^“*^ + *.. + a •, s t a 
O i II- n 

By a proper normalization procedure ["See Appendix 
equation (A-E-8) I , k , a and a can be made equal 
to unity. The normalized transfer function, say F(s) , 
is given by: 


1 

F(s) = 

s^ + 8^“^ + ... + a^_3_ s + 1 


(4.2.2) 


It is quite interesting to observe that the transfer 
function (4,2.2) can be modeled by means of a resistively 
terminated L-C ladder network f2B 1 1 The procedure to 
arrive at such a model will now be explained elaborately, 
following the references [25,26,27 * . 

Tile model sought is of the form as shown in 
Figure (4,2.1-a). Using the same notation as in the 
reference cEej , the transfer admittance function Y (s) 
and transfer function G-^^gis) for the model are written 
as: 


ha'®' ®a 
\(3) V yaa'®' 


(4. a. 3) 





T„(s) 




yig(s) 


G2+ygg(s) 
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(4.2.4) 


wliere G-g = 1/Rg and yj^gCs) and yggCs) are 
respectively the short circuit transfer and the short 
circuit driving -point admittances of the network shown 
by a rectangular block in Figure (4.3.1-a). If ^ 2 = 1 , 
then ^12 = ■^12- Then following Cauer’s method (Cauer’s 
first form) , the elements of the model can be determined 
quite easily. If 


^(s) = 


P(s) 


p^{s) + PgCs) 


(4.2.5) 


where 

P(s) = s^ + a^ s^”^ + ... + s + 1 (4.2.6) 

and p (s) and p (s) are even and odd parts of P(s) 

1 2 

respectively. Rewriting (4,2,5) as 

l/p^(s) 


F(s) = 

1 + P2(s)/p^(s) 

(4.2.7) 

and comparing (4. 

2.7) and (4.2,4) result in: 



= l/p^(s) 

(4.2.8) 


= P2{s)/p^(s) 

(4.2.9) 


y^gis) in (4.2.9) can then be expanded in a continued 
fraction to determine the values of the inductors and 
capacitors of the network in Figure (4.2, 1-a). The model 



He.T)rYt, FdTdds, Ohms 

id) 


network /NTERPRETATION for SCHWARZ CANONICAL FORM 
FOR SINGLEVARLABLE SYSTEMS 
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with all its elements known can now be represented, as in 
Figure (4.S.l-b) or as in Fig-ure (4.2,1-c) depending 
upon whether n is even or odd. 

Example : Consider 


F{s) 


T(S) 


From (4.2.9) 


yggls) 


Pg(s) 


1 

s^ + 2s^ + 3s^ + 4s + 1 

4 P 
S + 3s + 1 

o 3 

3s + 4s 


_ 1 


s + 


2s + 


s + 


2s 


Hence the proposed model for this transfer function is 
as in Figure (4.2,1-d)* 

As Routh’s stability criterion may be formulated 
through continued fraction '16, 28 ' it is quite reasonable 
to expect that the elements of the ladder network can be 
determined from the knowledge of the Routh array of P(s) 
in (4.2.6) ’.25J . Define H(s) as 


H(s) = 


a + a^ s^~® + 

1 o 


n n-2 

s + a^ s + . , 


(4.2.10) 
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H(s) 

His 

where 


From 

where 

P(s). 

D. 


can be expanded in a continued fraction. 

I = 1 


s + 


a 


1 s + 


b 




s+ * . 


(4.2.11) 


a 

b = ap - ~ 
1 2 a. 


‘l ~ ^3 " 


hi 


d = t 

1 3 c 


(4.2.12) 


. . .etc . 

[4.2.12) it is obvious that 

b^ = R 

1 31 

c = R (4.2.13) 

i 41 

. . .etc. 

R. .’s are the elements of the Routh array of 
3 

Define 
1 L, 




k odd 


k even 


(4.2.14) 


(k = 1 ,2 , . . . ,n) 
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wiiere L and C refer to tile inductance and capacitance 
respectively. Then from (4.2.11) and (4.2.14), it is seen 


that 

R 

, (k = l,2,...,n) (4.2.15) 

^(n-k+2)i 

Thus the first column of the Routh array can be 
effectively used to determine the elements of the ladder 
network. If the system is stable, all the (k=l ,S, . . . ,n) 

will he positive. In other words, stability is a necessary 
condition for obtaining a network f Figure (4.2.1-a)J which 
is physically realizable. 

From (2.3.10) and (4.2.15), it is very simple to 
see that 


b 

n 


1 


D. D , 
1 i+1 



(i = 1,2 


n-1) 


(4.2.16) 


It will be next shown that a proper choice of 
state variables will yield the system matrix as any one 
of the three Schwarz canonical matrices (S^ , Sg and S^) 
for the system which is modeled by a network as in 
Figure (4.^.1). Although in the followir^ treatment n 
is assumed to be even, an exactly similar treatment; is 
possible for ’n odd’ case also. 

Consider Figure (4.2.1-b). If charges associated 
with capacitors and flux-linkages associated w’ith inductors 
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are chosen as state variables, the system equations 
can be written as: 


m 

s 

where 

s 


D s + e V 



3 ’ ^ 4 ’* 


} 


0 - 1/0 0 

dj 

l/L^ 0 -1/Lg ... 

0 1/Gg 0 


0 

0 

0 



0 


(4 


0 


0 0 


and e = (1, 0, ... , 0) 


T 




Choose the 

state vector 

w such 

that 


- >1 





- ( ^ 1 

^ k=l 

> 

)"7i 

i odd 

W 

X , 






^ k=l 

/ 

a, 

1 

i even 


\ 

(i 

= 2,3,. 

. . , n) 


In other words, seek a transformation 
w = T s 


r 


(4 




(4 


where 


T 


diag. 



-h°8- 





(4 


.2.17) 


.2.1b) 

, 2 , 1 ? * 

. 2 , 20 ; 
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It may bo easily verified tliat 

4 . 


and 


T D T' 


To = d 


es 


where 


(- 1 , 0 , 0 ) 




(4.2.21) 


Therefore 


w * Sg w + d V 


(4.2.22) 


Thus when the state variables (i=l ,2, . . , ,n) are 
chosen as in (4,2,18), the system equations will' be 
given by (4,2.22). 

Consider now the state variables (i=l ,2, , ,, ,n) 


given by 
/ 


‘li 


yi=( 


(5 i 

Yi 


i even 


i odd 


/ i- 

^k=l ^ 




(i=l ,2, , , . ,n) 


or 




(4.2.23) 




JJ 3 


(4.2.24) 


where 



? 

With 

where 

that 

n-i+l 

i .e. 

where 

V = 

1 


= dlag. ; (- h), ( 


54 




•), ( 


h <^8 1-3 ’’"7 


(4.2.24), (4.2.17) becomes 
y = B y + g V 


(4.2.25) 


(4.2.26) 


f ~ ( - ~ 5 0 , . . .,0 )^ 

1 


(4.2.27) 

Pick the state yarlables (1=1 .S , . . . ,n) such 


/I 

k=l 


D, 


7 


^ i- 

k=l 

(■^~^ 3 2 > • • . , jtl) 


"7 


I 3 i odd 


, i even 


(4.2,28) 




r = V s 


(4„2.2C) 


0 

1 


(L2Cglg. . . t) 


Jl> 


L C 
1 2 


(4.2.30) 



mzix (4.ci.kigy 


where 


( 4 . 2\ 1 7 ) is transformed to 
r = Sg r .+ h V 

h = ( 0 , 0 , ... , 0 ,- 1 /L;L^'^ ' 


(4.2.31) 

(4.3.32) 


Can this investigation help in getting the 
transformations between and Sg and between 
and §2 ? Let these transformation matrices be 
respectively and M 3 . Clearly 

■^ = Us = U w = Mg w 

JL • 6 • 

Mg = ? (4.2.33) 


u and T are diagonal, Mg becomes diagonal. It 
is readily seen that 



(4.2.34) 


If it is reqioired to know (TM^)^^'s (i=l ,2, . , . ,n) in 
terms of then equations (4.2,16) may be solved 


to get: 
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n-i-l 

fn=0 ,2 ,4, . . . 


n-i 

\ •u 

“^n-k 

k=0 5 1 ,2 j ^ , 




n-i-l 


/ 


I. 


(b )' 
n-m 


m=l ,3,5,... 


n-i 


-uZ ^n-k 
k=0 , 1 , 2 ,... 

(i=l,2, .. . ,n-l) 


D = 1/b^ 

n ' n 


(n-i) odd 


(n-i) even 


(4.2.35) 


Using (4.2.34) and (4.2.35), (Mg)^^ (i=l , 2 , . . . ,n) can 
be computed from b^’s (i=l ,2, . ,n) . The diagonal 
transformation Mg is simpler than the one (Tg) obtained 
in chapter II See (2.2.16) . As 

y = Us = UV“1 r , (4.2.36) 


Mg can be detennined from 

= y (4.2.37) 

It is quite simple to see that Mg is a cross-diagonal 


matrix and that 
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(Mg) 


i (n-i+1 ) 


0 when j ^ n-i+1 

(- 1 ) 1*1 

{i=l,S, .. . ,n) 


{4.8,38} 


J 

Thus it has been shown that each of the mathematical 
models (4.3.28), (4.2.26) and (4.2,31) can be thought of 
representing the state equations of a ladder network in 
Figure (4.2.1) , when the state variables are appropriately 
chosen. In fact, each state variable is proportional to 
a corresponding variable (charge or flux linkage) 
associated with an element (capacitor or inductor) of the 
ladder network, the proportionality constant being related 
to the parameters of the network. 

It may be observed that the mathematical model 

(4.2.22) or (4.2.26) or (4.2,31) has been 'realized’ by 
a physical network shown in Figure (4.2,1). in this 
sense, it is a synthesis problem. It may be expected 
that such an interpretation may facilitate analysis too. 
Consider any physical system which can be modeled by a 
ladder network in Figure (4.2.1), The state equations 
for such a system can be directly written in the form 

(4.2.22) or (4.2.26) or (4,2,31). If the system 
transfer function is given as in (4.2,1), then by 
employing the procedure outlined above, the state equations 
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Suppose tile given system is described by 

X = A z + b u (4.2,39) 

where A is any arbitrary n x n matrix. Reference f 21^ 

has given a transformation (See also chapter III) which 
takes the system in (4.2.39) to the one given by: 

y = B y + f u (4.2.40) 

where 

f = ( 0,0,. ..,0,1 )T (4.2,41) 

provided (4,2.39) refers to a controllable system. Equation 
(4.2.40) is same as the equation (4.2.26) except for the 
input vector f. As f is different from g , the network 
corresponding to (4.2,40) is not the one in Figxire (4.2,1) 
and it is given in Figure (4.2,2). 

4.3 NETWORK INTERPRETATION FOR 'TILTITARIaBLE SYSTE^'IS 
IN SCHWARZ FORM 

Schwarz canonical form for multivariable systems 
was introduced in chapter III, It has been shown that 
any controllable mtilti variable system 

X = ^ 5 + F y (4,3,1) 

can be transformed through 

y = P X (4.3.2) 

to a system 

■ ‘ 

y = B y + F V (4.3.3) 
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The system In (4.3.3) consists of coupled subsystems each 
of which, being in Schwarz* form. Thus, each of these 
canonical subsystems can be associated with a resistively 
terminated L-C ladder network with voltage or current 
drivers at appropriate locations in the ladder network. 

To account for couplings between canonical subsystems > 
dependent drivers will appear in the networks. These 
ideas may be better explained with respect to an example. 


Example ; Edr convenience, consider a multivariable 
system (4.3.1) which is already in phase -variable form 


[i.e 

. A 

It 

1 

and 

P = E. 

See 

seotion 


with 





\ 




- 



0 

1 

0 

1 0 

1 

0 

0 

\ 

0 

0 


0 

0 

1 

I 0 

0 

0 


0 

0 

A = 

-2 

-3 

-1 

1 -4 

_ 1_ _ _ 

-2 

-1 

and E = 

1 

- 5 


0 

0 

0 

1 0 

~r 

0 


0 

0 


0 

0 

0 

1 0 

0 

1 


0 

0 


i -1 

-2 

-3 

1 

-9 

-2 


0 

1 





1 



i 




By the method discussed in chapter III, the transformation 
matrix P is found to be 


1 

0 

0 * 

1 

0 

0 

0 

0 

1 

0 1 

0 

0 

0 

2 

0 

1 1 

0 

0 

0 


. — - 


_ _ 

- - 

- 

0 

0 

0 1 

1 

0 

0 

0 

0 

0 j 

0 

1 

0 

0 

0 

1 

0 

3 

0 

1 
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Tiie input transformation matrix G- See (3.5,4) can 
be written as 

1 -5~ 

G = 

_ 0 1 

TJie transformed system becomes 

y = B y + F V . (4.3.3) 


Where 



Note that the coupling terms in a and B are different 
( owing to the transformation) . Write down B ih the 
following form: 



whore the values of and g^’s are as given 

above "^See (4.3,3) J , The ladder networks associated 



85 ^ 11 ^ 12^13 


^^ 11 *^ 12 ^ 13 ^ 


*4^11^12^13) T f*6^11^12^13 


21 '^Z 2"23 


Dependent drivers due to coupling terms ✓ 

, 84» 85 86 ( *2 4 1^2 2 ^2 3 

23 ^^2i*-2Z^23) V2 L,,C,. 


g 11 ^ 21 ^ 22^3 


* 3 ^ 21 '"2 2 ^2 3 


'11 '^ 12^13 


Dependent drivers due to coupling terms 

gj. §2 g 3 


NETWORK INTERPRETATION FOR SCHWARZ FORM 
FOR MULTIVARIABLE SYSTEMS 


FIGURE (4,3,11 
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with the subsystems in (4.3,3) are given in figure (4.3,1), 
The various quantities ^13 ^21 j 

Cgg and Lgg can be computed respectively from ^12 

and b^g and tgg ^ud bg^ using the relations in 

(4.2.16). 

4.4 COWCLTBION 

Some possible uses for the network interpretation 
introduced in this chapter have been given already in 
section 4.2. Recently Israel Navot L293 has considered 
the synthesis of certain subclasses of tridiagonal matrices 
with prescribed eigenvalues. For stable eigenvalues 
(with negative real parts), Schwarz matrix can be quite 
simply (and uniquely) synthesized. The network 
interpretation given in this chapter can then be readily 
used to "realize" physically the synthesized Schwarz 
matrix. This results in resist! vely terminated L-C 
ladder network. If synthesis in terms of R-L-C ladder 
networks is contemplated, the necessary methods are given 
in [^27,30, 31J , The network interpretation discussed 
above will be used to justify the procedure employed for 
simplifying large dynamic systems in chapter V. 





CHAPTER - Y 

A IffiTHOD FOR SIMPLIFYING LARGE DYNAMIC SYSTEMS 
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5.1 INTRODUCTION 

One of the promising features of state-space 
techniques is that they are applicable to systems (single- 
variable or multivariable) of any order. However, exact 
analysis of large dynamic systems (a chemical process 
control system, for example) demands excessive computational 
time and cost. Thus it is usual in control systems 
practice to analyse high order systems through approximate, 
low order systems. For instance, aircraft systems are 
often analysed through second order or third order models. 

The computational aspect may now be considered. 

In a typical computational routine, the integration 
interval is to be of the order of, or less than the 
smallest time constant of the system. Thus, for a system 
having some eigenvalues much larger than others (larger 
eigenvalues mean smaller time constants), such a routine 
calls for short time increments which would mean an 
increased cost of computation as the number of steps in 
the computation will be large with short time increments, 
in such cases, the modes associated with ’’large” eigenvalues 
(high frequency modes) may be neglected, thus reducing the 
computational complexity. In this work, such an 
approximation will bo called as dominant-mode approximation, 
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even though, it is customary to define a pole (or a pair 
of poles) to he dominant if it (or they) lies (or lie) 
nearest the imaginary axis of the s-plane L32,39j/ . 

Recently there have been many papers dealing with 
the simplification of large dynamic systems _33,34,35 ,36 ,37 
38 J, , Davison’s method l33] requires the knowledge of the 
eigenvalues and eigenvectors of the original system in 
order to be able to determine the system and input matrices 
for the simplified system. Gustafson [SSJ employs 
truncated and associated transfer functions to approximate 
any order system by a second order system. These methods 
do not require the computation of the eigenvalues and 
eigenvectors of the original system. The matrix 
generalization of the truncated transfer function method 
through phase-variable canonical form is given in [36 J , 
Here, the Schwarz canonical form is used for the matrix 
generalization of the associated transfer function method. 

All simplification procedures which retain the 
dominant eigenvalues (with or without modification to 
account for the neglected non-dominant eigenvalues) suffer 
from the serious drawback that the approximation made is 
not good for all control inputs. This can be best 
illustrated by means of a simple example. Consider a 
system described by: 
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The eigenTSlues of this system are -i and -100* 

Wxon uCt) is a step function, 

X (t) = 1- 0.99 + — Q-lOOt 

99 

4. dominant-pole approximation of the above system l40j 
yields 

z - X + u(t) 

where * denotes that the variable is only approximate, 
Ihat is, the mode corresponding to -100 has been 
neglected. With the same unit step input, 

x*(t) = 1- e”^. 

Mote from the expressions for x!^(t) and- x (t). that 

** 3 ^ 1 , 

X2(t) - and X (t) are reasonably close., Gonsideribg 
the s-plane, the step input u(t) introduces a pole 
at the origin close to. -1 and hence the approximation 
is reasonable. Suppose u(t) = This input 

function introduces a pole at -99 close to -100. Now, 
with this u(t ) 

' (t) . —122 ,-t . 100 g-99t wo e-wot 

1 98 z 99 98 99 

■■ v 

4 ' 



and 


1 

98 
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x* (t) 


98 


The strength of the mode associated with -100 in the 
expression for 2 :^(t) is now much larger than that 
associated with -1. However, the former mode decays more 
rapidly than the latter one. Thus, for t» 1/100 second, 
X (t) is "close” to x^(t). If the transient response 
in the interval {/o, O.OOIJ is of interest, then the mode 
of -100 becomes dominant. The simplified tion procediice , 
to be described in this work, also suffers from the above 
drawback. However, after "sufficiently large” time, the 
approximate transient response would hopefully approach 
the exact transient response. 

Section 2 of this chapter states clearly the 
type of systems considered for simplification and outlines 
the simplification procediire. Justification of the 
proposed method is given in section 5. Sections 4 and 5 
illustrate the method through several examples. Section 5 
also discusses some difficulties in reducing systems 
having complex eigenvalues. Problems for future 
Investigation are posed in section 6. 


5 . 2 SIIIPLIPICATIO N PROCEDURE 

The system (with all eigenvalues in the left-half 
plane) to be considered for simplification is described 
by the transfer function 
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G(s) = 2jL£i = 1 ( 5 .S. 1 ) 

u(s) + ags^”^ 

^ few remarks on tills transfer function are now In order. 

G(s) in (5,2.1) is actually a normalized transfer function- 
In tile usual form (not normalized), it can have in the 
numerator some gain constant k and the coefficient of 
in the denominator need not "be unity. In such a 
case, the transfer function can te normalized to the form 
in (5.2,1) through a simple normalization procedure [ See 
Appendix II, equation (A.2,5)J , Next, G(s) does not 
have zeros (no numerator dynamics). This constraint can 
be quite easily relaxed. Tor details refer to Appendix III. 
Finally, throu^ a simple choice of state-variables, the 
system in (5,2.1) can be put in matrix form. Define 

x^^ = X , Xg = X , (5,2,2) 

The equations (5.2.1) and (5,2,2) yield 

01 0 ...oo^jxj^- fo 

0 0 1 ... 0 0 Xg 0 

0 0 0 ... 0 0 x„ 0 

• o 

• • • ••••• • * 4 ^ • 

9 0 0 ••••• • • 

00 0 ... 0 1 X^_^ 0 

^n~®n-l “^n-2 *** “^2"^1 \ 


(5.2.3) ^ 




• • » 



or 


sc 

X = A X + f U (5 .2.4) 

where x and f are n-vcotors and .i. is an n x n 
matrix. Thus the system, irsteaa of being given in the 
form (5.2.1), could as well have been gj^^-en directly in 
the form (5.2.4). In fact, it need noc be in phase- 
variable form but can be in any arbitrary form. In 
such a case, a suitable co-ordinate cransf ormation [_15,22 j 
is to be applied to (5,2.4) so that the transformed 
system is in phase-variable form'^. laving pointed ouc 
this, hereafter a will be as=!Uiiod to be in phase-variable 
form. 

The simplified system is given by 

X* = A* X* + f^ u (5.2.5) 

Where x* and f* are -^.-vectors ( i < n) apd en 

matrix. Further, f*' has the first (^^'-l) rntrios 
ecual to zero and the last one :s aifto-^ont from zero and 
ii* is in phase-variable form. 

The system in (5.2.4) js tr a nsf orricd to its Schw-rz 
form, simplification is perfo’^med at -chis stage and thG 
simplified system is then transformed back to the form 
in (5.2.5)., This is the philosophy of the simplj ficotion^ 
procedure. This will be now olaboj''tely explained. 

*Such a transformation requires that the system mist bo 
completely controllable., 



Seek a transformation 


z = P X (5.2.6) 

■which takes (5.8,4) to the form 

i “ ? 5 ^ (5,8.7) 

where B is in Schwarz form i.e. 



0 

] 

0 ... 

0 

0 

* 



-bi 

0 

1 ... 

0 

0 




0 

-bg 

0 ... 

0 

0 



B = 

« 

• 

. ... 

• 

« 


(5.2.8 


« 

0 

* 

0 

. ... 

0 ... 

0 

m 

-b 



i 




n-1 

n 



Compare now 

b - 
n-1 

with 

^n-2* 

Is 

b -r 

n-1 

much larger 


than o The comparison is said to be successful 

■u n-fc- 


if ^ 10. This constant 10 has been found to 

^ n 2 

be ~ suitable for many problems. V/hether this 

comparison is successful or not, go to the next comparison 

i.e. compare ^n-3* i-th comparison 

is successful if 77^^^ ^10. If during the process, 

b . ^n-i-] 

— becomes less than 1, stop comparing any further. 

Td 

n-i-l xf this situation does not happen, comparisons 
are to be carried on till the step in which ^>2 aud b^ 
are compared, Now ;he statement of the result is in 
order. If no comparison turns out successfully, infer 
that no simplification is possible. If i-th comparison 
is the last successful one (some or all comparisons before 
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this i-th successful comparison could have been unsuccessful) , 
then the system in (5.3.7) can be simplified and the 
simplified system will have the size eCLual to (n-i) . Let 
^ = (n-i). The simplified system in Schwarz form is 
given by 


B* z* + f* u 


where 


= 


-b* 

1 

0 


0 

0 


1 

0 


0 

0 


0 

1 

0 

• 

0 

0 


0 

0 

0 


0 

0 

0 


0 1 

■t-1 



if 

O 

o 

b* 

= , 

1 


bh 


— '•"**' " 

H*! 

■n 

_ 1 




b/ 


+3 


+3 * * • 




'*^-^ + 1 ^-^43 •” 

Transform (5.3.9) to (5.3.5) through 


(5.2.9) 


(5,2.10) 


(5.2.11) 

(5.2.12) 


(5.2.13) 


z* = P* X* 


(5.2.14) 
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Several comments on this method will now be made. 
First, P and P* can be very easily formed tiiroiJgJi 
the procedure given in chapter II. Second, the above 
method is a matrix generalization of the associated 
transfer function method discussed in [bsJ . Reference 
i_25j considers a transfer function G(s) given in 
(5.2,1) and reduces this to: 

G'(s) = (5.2.15) 


where “the elements of the Routh array of 

the denominator polynomial of G(s). To illustrate 
the idea that the proposed method is a generalization of 
the associated transfer function approach, let us write 
'down the simplified transfer function (of order ^ ). 


G*(s) = 

u(s) 


or 


G’^(S) 


I 




-J- E 


'f-l 


+ • • • +R 


(n+l)l 



1/R 




R 




/-I 






(y)-l+0 I 

(5.2.16) 


(5,2,16) can bo recast as 



f u 


Mb) 


f 


u 


(s) 


wJierG 

Q* 


= 


t *<* a, 

s'+a^s +***+af 


i. 


E 

"e 




E 


(n+l)l 


E 




y 


(5.2.17) 


(5.2.18) 


Note here that the factor miilti plying u(s) in (5.2,17) 
l^i.e. is equal to p in (5.2,13). The 

relations (5. 2, IS) may be easily obtained from (5.2,18). 

Thirdly, we shall illustrate as to what such 
an approximation means. For convenience, consider a 
fourth order system given by 


u(t) (5.2.19) 


• 


■n. 



- 

^1 


0 10 0 


^1 


0 

^2 


-b^ 0 1 0 


^2 


0 








^3 


0 -bg 0 1 


"3 


0 

• 1 

1 


0 0 -b -b^ 


z . 


1 

■ ^ J 


5 4 


4 


— ^ 


Suppose (bg / bg]^10. So it is inferred that this 

system can bo approximated by a system which is at least 

one order less than the system in (5.2.19). Instead of 

comparing b with b right now, the system in (5.2,19) 

1 

may be reduced by one order and then ^ 




computed. As per equations (5. 2. '9) - (5.2.13), tJbie 
simplified system is given by 
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(5,2,20) 


This could have been obtained if we had equated z in 
(5,2,19) to zero i.e. 

= 0 = Zg -b^ ^4 ^ 

24 = -(^3 ^ 4)23 +(1A4)u (5.2,21) 

Substituting (5.2,21) in (5.2.19) resiilts in (5.2,20), 

Next, if bg /b^ ^ 10, then a similar procedure can be 
used to get the 2x2 simplified system. This viewpoint 
helps us to approximately estimate the variable 
corresponding to the neglected mode (in this case, z^) , 

As we have forced z^ to 0 , the value of found 

through (5.2,21) may not represent a good approximation 
of actual z^. However, for the sake of completeness, 
let us investigate the following. Suppose the output 
equation of the system in (5.2,4) is given as 

y = C X + d u (5.2,22) 


A 

Denote y as the approximate estimate of the output y , 
to bo evaluated from the simplified system in (5.2,5). 
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Denote z as tJie approximate estimate of z in (5.8.7). 
Tlien, it may be verified that 

(5. 2. S3) 



2 

= 5 2* + e u 

where 



5 nz<' 

=r 

-In. -^2 



‘■nxn jiur 


== 

0 , 1 ^ j 

(K, ) , . 

“1 11 

= 

1 > i = 1,2,.. .,ii 

(Kn ) 

^ (^+1)(^+1) 

(k ) 

-1 2 ) 


H ^(l+ 2 )*». 


^ (-^-+1) "b ( t+ 3 ) . . . 

= 


(Ui)(i^i) 

= 

1 ^ {^+i-2)(y+i-2) 



i = 3 , 4 ,...,n-y 

is 


|^(n-£+l) ( ) ^(n-y/+l)l 

1 

— ^ 

( 1, 1, ..., 1, 1 )!' 

-nxl 

"ij 

r= 

^Inxl 


0 i 3 

hi 

= 

1 , i = 1,2,...,^ 




JS 


(^+1) (-(i+1) 


^(^+l)^(^+3) **• 


H 


(^+i)(-^+i) 


§1 




(5.2.S4) 


D > i even 

^(■^+1-8) (t+i-2) » i o<3d 
i = 2,3, ... ,n- t 


0 


'xl 


“(n--? )xl 


- (n-^ )xlj - 

( 1 , ... , 1,1 )^ 
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Again, as mentioned above, , z , . . . , 'z 

will be ’’very” approximate, as we bad equated the 
derivatives of these variables to zero. Now, as 

X = p-1 z , (5.2.25) 

X = ?-l z = p-1 K z* + P“1 e u (5.2.26) 

(Note carefully that x and z are n-vectors and x* 
and 2 * are -vectors). Xfeing (5.2,22) and (5,2.14) 
y - C* X* + d* u 

where 

[ 

C* = C P-1 K P* ^ (5.2.27) 

and 

d* = C P"^ e + d 

It is worthwhile to repeat that ^ will be ’’very" 
approximate and will represent a reasonable approximation 
of y Tactual value determined by (5.2^22)] if C in 
(5.2.22) is such that y depends only on the first 

i 

*' -components of the state vector x. 

Fourth, we shall make use of network interpretation 
concept discussed in chapter IV to investigate this 
approximation process. Again, for convenience of 
explanation, consider a fourth order system given in 
(5.2„19). This system can be modeled by the ladder 
network shown in Figure (5.2,1-a). The state-variables 
associated with the network (marked in the Figure) and 







78 



as; 



... (5.S.29) 

Suppose (bg /bgj ^10. As bg = 1 /( 130 ^) and 

bg = l/iCgLg) , this inequality implies C^/Cg 10 . 

Now, equate q^ to zero, thus find i-ii terms of 
y 3 and u and substitute back in (5.S.19) to get 



(5.E.30) 
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Equations (5.2.30) correspond to a ladder network 
shown in Eigure {5.2.1-b), The Schwarz form for the 
network in Eigure (5.2,1-b) can be readily written as 



and this is exactly the same as the one obtained through 
(5.2.9) - (5.2.13) with and n=4. Remember that 

q , above, was equated to zero. Thus a rough estimate 

41 : 

of q^ can be found from this equation. It may be 

verified that this woiiLd, after applying the transformation 

2 = ^ ^ — j^Seo equation (5.2.28)J , result in 

^ 1 2 '^ 3^4 

(5.2.21). 

A careful observation of Eigures (5 ,2,1-a) ^and 
(5.2.1-b) reveals that the network in Eigure (5,2,1-b) 
may be obtained by deleting the component from the 

network in Eigure (5.2.1-a) and then transforming the 
current source in Eigure (5.2,1-a) to a voltage source 
as in Eigure (5.2.1-b) . This viewpoint [ that is, 

neglecting the (n--^) dynamic components near the resistor 
end of the ladder network when the simplified system 
size is determined to be tj will be heuristioally justified 
in the next section. 

Fifth, it is to be emphasized that the proposed 
method does no-^: require the computation of eigenvalues 
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major advantage over Davison^s method i_33j. However, 

the method proposed here shows, at present, some difficulties 

with systems having both real and complex eigenvalues* 

These are elaborately discussed in section 6, It is hoped 
that in future, these difficulties would be resolved. 

Finally, a couple of more favourable comments will 
be made. Construction of Liapunov function for the reduced 
or simplified system is q.\iite trivial. Denoting V* as 
the required Liapunov function, 

Y* = x*"^ ( L* P* ) X* (5.2.31) 


where 


!» . ( 5 . 2 . 3 £) 


Davison’s method results in a reduced system matrix A* 
which retains all the dominant eigenvalues of A without 
modification. But the method proposed here modifies the 
eigenvalues in order to account for the neglected modes. 
It is believed that this method would yield improved 
approximation. Also note that if the original system is 
asymptotically stable*, the simplified system is also 
asymptotically stable as all the bf's (1=1 ,2, , . . 
are 0, 


* In this Work, only asymptotically stable systems 

are of interest. This fact has been mentioned in the 
beginning of this section. 
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5.3 JUSTIFICATION OF THE SBIPLIFICATION PROCEDURE 

TMs section alms to show henristlcally that 

whene-ver the given (original) system has 1 eigenvalues 

which are ’’much larger*** than the smallest eigenvalue 

of the system, the i-th comparison will he the last 

successful one (refer to section 5,3). As the h * * s 

i 

(i=l ,3, . . . ,n) of the Schwarz matrix depend nonlinearly 
(through Hurwitz determinants) upon the coefficients 
of the characteristic polynomial and these latter 
coefficients themselves are related to eigenvalues in 
a complicated manner (Any book on theory of equations 
will give these relations), a rigorous justification 
of the above statement may not be possible. Instead, 
a heuristic reasoning will be provided in what follows. 
Generalization, based on the experience gained by 
considering low order systems, of the results to any 
n-th order system will be made. Before proceeding 
with the heuristic justification, it is to be 
mentioned once again that the whole simplification 
procedure works very well with real eigenvalue syjstems 
(repeated or distinct) , but encounters some difficulties 
with systems having real and complex eigenvalues* 

*An eigenvalue is said to be large if its magnitude 

is more than say k times that of the smallest eigenvalue 

of the system, where the value of k has to be decided 
depending upon the accuracy desired. A convenient fiauve 
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Consider a 3x3 system. It is assimaed that the 

eigenvalues /■, and ^ of this system are real and 
1 2 o 

stable, and are such that one of them (magnitude only), 
say is much larger than the others i.e,/\g and 

This implies that if 



^2 

/ 

(5.3,1) 



' i 


then, 

k^ and 

kg « 1. Let r^ = - (j 

.=1,2,3)*. If 

the characteristic polynomial of this system is given hs 


P(s) 

rz 2 

= S -j* Q, B + ^oS 3|-? « 

1 3 ' 

(5.3,2) 

then 


= ag/a^ ; hg = ag 




and b„ = 

3 i. 

(5.3.3) 


From the relations between eigenvalues and the coefficients 


of P(s), (5.3.3) can be rewritten as; 

^ 2 
b^ = kg r^ 


(k^ + 


= (1 + k^ + k2)r^ 


( 5 . 3 .. 4 .) 


*The relationships between the ’s and 
the coefficients of the characteristic ^ polynomial 
Involve powers of (-1). So to avoid such factors, 
are defined and now the resulting relations will be 
simpler [^4:1 J . 
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As 


and 


kg << 1 , it is obvious from {5.3.4) tiiat 



»1 


(5.3.5) 


thus resulting in one successful comparison. WJien 

~ ^2 (i-e. )\2 ^Ad are equal and they are 

much, less than /\^), (5.3,5) still holds. Suppose k^ 
is comparable to 1. Then for kg 1 


^2 

ki kg 





» 1 


Therefore, even in this case (5,3.5) holds. 

Let us now consider a fourth order system. With 
the same notations as used in the previous case, the 
coefficients of the characteristic polynomial can be 
expressed in terms of the r^^’s (i^l ,2,3, 4)» as* under: 


a 


1 


a 


2 


a 


3 


a 


4 


Also, 


b 

b 


1 

2 


b 


3 


b 


4 


[l + (k^H- kg + kg)! r^ 

|^(k^ + kg + kg) + (k^kg+ ^2^3 ^ ^2^3 li 

[(k^kg ^ ^ 

(k^kgk^) r^ 


34/(82 -^3/^1) 

^3/^1 “■^ 4/^®2 -^3/^^ 
^2 


(5.3.7) 
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Now 

^3 "" Uk3_+k2+k3)+(k]_k3+k]_k3+k2k3) 

”*> 

-(k^kg+k^kg+kgkg+k^kgkg) /(l+kj^+kg+k^lJ 

... (5 ,5.8) 

^2 "" “ ^4 “®3'^®4^ 

r" 

2 *' 

= ; (k^kg + k^k^ i- kgk^ + k^kgk2)/( l+k^+kg+k^) 

k^k2kg,(l+ k^ + kg + k^) 

(k-j^+kg+kg) ( l+k2^+’kg+k3+k-j^k2+k2k2+kgk3*' k^kgk^)^' 

...(5.3.9) 

Having written down tiie exact expressions for tg and bg, 
let us investigate the following cases. 

Case I : Suppose k^, kg and k^ << 1. That is, As* jig 
and are much smaller than /\^. Then from (5,3.8), 

bg = r^ ( k^ + kg + kg ) (5.3.10) 

In arriving at (5.3.10) from (5.3.8), two approximations 
have been made. They are; ■ * • 

(i) k^kgkg is neglected in comparison with 

(kl+k^-k^) and * kj V W : 

(ii) (k^kg + k^kg + kgkg) is neglected in comparison 

with (k + k + k ) : 

12 3 

These approximations are reasonable as k^, kg and 
kg are :'C 1. Neglecting similarly k^ kg kg and 



(k-j^kg + ^1^3 ^2^3^ comparison with (l+k^-** kg + kg) 

in (5,3.9) wo get 
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b = r^ 




h ^ 'T J^. ^ -TiY 2 Til. g (k^-^-kg+kg) (l+k^+kg+kg) 1 


l-fk^+kg-i-k 
(5,3,11) can be rewritten as 


(5.5.11) 


r^ 

■^1 

( ■ ^1^3^3 / 


{ . c,. +k,k„+k_k„ "• / 

L " 1 3 S 3 (kj+kg+kg) j 


(5.3.12) 


Thus , 


^2 < 


l+k^+kg vkg 


{k.k„ + k.k^ + k„k„) 
12 i 0 20 


(5.3,13) 


From (5.3,10) and (5.3,13), it is easily seen that 


bg (kj^+kg+kg) ( l+k^+kg+kg) (kj^+kg+kg) 

J: > — — -_> -- — (5.3.14) 

8 

bg (k^kg‘>'k^kg«5'kgkg) (k]_kgfk2kg+k2kg) 


AS K J k_ and k^ are 1 , (5 ,3.14), implies 

1 O t.> 


^3 


(5.3.15) 


Recapitulating, we had a fourth order system with one 
eigenvalue much larger than others and hence we obtained* 
a successful comparison given in (5,3.15), 
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Case II : Suppose k^, kg and k^ are all not muoii 

less than one i.e. for instance k]_ = 1, but kg and 

k_ may be much, less than one. Even In this, going 
5 

through the steps given in Case I, it is easily seen 
that bg is larger than b^ (bg>bg) but not necessarily 
much larger than bg. Also, it is readily seen from 
(5.3.7) that 

1 + bg / b^ — ^^ 3(^2 “^■ 3 /^ 1 ) (5.3,16) 

Substituting for a-»s in (5,3,16) and rearranging yield 

T ^2 j, 1 1 . 1 ' I ^^l"^ ^2 “^^3^ 

1 4 — » = '' 1 + + 4- I I 

^2 l(l+k^ H-kg 

[i+kl-<-kg+kj+kilCg+k^kg+kj.kj-k.':gkj(l+ki+lC3+kg)] 

... (5.3.17) 

Now making use of the fact that ~ ^ ’ ^2 ^5 ^ 

and after some simple manipulations we obtain 

1 + bg/b^ = (kg + kg) / 2k2k2 (5.3.18j 

AS kg and k^ 1, (5,3,18) implies 

bg/bi >) 1 (5,3,1c) 

(For example, if kg = 0.01, k^ = 0.02 and k 2 _ = 1, 

bg/b^ becomes approximately 75). Summarizing, we had a 
fourth order system with two eigenvalues comparable to 
each other in magnitude and much larger than the other 
two and this resulted in a situation wherein the second 
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comparison was tiie last successful one. 

Case III ; Suppose the system has two distinct real 
eigenvalues and one repeated eigenvalue of multiplicity 2. 

Two possibilities are now open to us. In one, the repeated 
eigenvalues have to be retained in the reduced system. In 
the second case, the repeated eigenvalues have to be neglected, 
From the results of Case I and Case II, it may be easily 
seen that in the first possibility, two successfxil 
comparisons will be obtained and in the latter, two 
successful comparisons need not occur but the second 
comparison (bg/b^>' 1) will be successful. Thus in either 
case the proposed method will reduce the system by two order. 

In a general n-th order case also, the same 
procedure can be applied to verify the validity of the 
proposed method (of course in real eigenvalue case only). 

As expected, the expressions for bj_’s now become quite 
involved, and thus are not given in this section. 

In section 5.2, it is asserted that whenever we 
have the i-th comparison as the last successful one, the 
i dynamic components near the resistor end of the ladder 
network can be neglected. Let us heuristically justify 
this assertion. In view of the discussion presented above, 
we have thus i eigenvalues which are "large" compared to 
others. 
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Gustafson 44 _! states tJiat the "high frequency 
modes" accumulate near the resistor end of the ladder 
network (See section 5.2 for the Figures) and thus 
noglecting the dynamic components near the resistor end 
only eliminates the high frequency modes. (Note that the 
method proposed in this work gives a way to determine 
how many such dynamic components near the resistor end 
are to be neglected). Let us in what follows attempt 
to "prove"* Gustafson’s statement at least in real 
eigenvalue case"’’. 

For convenience of explanation, consider a fourth 
order system. The coirresponding ladder network for this 
system is shown in Figure (5,2,1-c). Now consider the 
Figure (5.2.1-d). Looking into the port 'ah’ from right, 
we have L^ and this resonates at the frequency = 0. 
Look pow into the port 'cd'. We have L^ and Cg in 
parallel. This parallel combination resonates at a 
frequency, say, fg and fg > 0. Let us now look into 
the port ’ef’. L^ and Cg are in parallel and this 

*This heuristic proof was suggested to the author 
by Professor S, Fenkateswaran, Department of Electrical 
Engineering, Indian Institute of Technology, KanpUr, 

. "''it is to be observed that Gustafson’s statement 

is for any general system whether it has real or complex 
eigenvalues. 
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coffl.'bination is in series -witli L^. At frequencies higher 
than fgj the parallel combination of and Cg becomes 

capacitive and this capacitance (call this as residual 
capacitance Cp„) is less than Cp L, neutralities 

^ vv\A. ^ 

a part of Og. In fact, at f^, the effective iepedeice 
of in parallel with Cg is equal to zero). It may 

be observed that for systems with real eigenvalues 
/ ii . > ( i=l , 3 , . , . ) and ,4 , , , , ) , 

Thus and Hence the resonating 

freq.uency of the circuit located to the left of 'ef*, 
say fg, is larger than fg l.e, fg ^ fg. Similarly 
we can show that the resonating fretquenoy f^ of the 
circuit located to the left of ’gh* is larger than fg. 
Hence we have 

f > f \ f \ f (5.3.20) 

4 ^ 3 ^ 2 ^ 1 


These frequencies roughly correspond to the eigenvalues 
of the system. Thus we infer that the dynamic elements 
responsible for the high frequency modes are located 
near the resistor end. The arguments presented above 
can be very easily extended to ary n-th order system 
with real eigenvalues. 

Is ’Gustafson’s statement true for systems with 
both complex and real eigenvalues ? For such systems, 
the statement that *^i ^ 

necessarily true. Let us give an example to illustrate this. 
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Consider 


G(s) 


+ o.ls^ + 100. Ols^ + 10.1s + 1 


For this system, = O.IH, = lOH, Cg = O.IF 

and = lOF. Thus in this system and Cg^C^, 

Suppose in the case of the fourth order system 
i.e. It is said above that can 

be neglected. Are we ;iustified in doing so ? The answer 
is yes if we can show under the condition C^<(<Cg that 
the impedance looking into ’ef’ is much less than the 
impedance due to 0/ [See Figure (5.2.1-d)3 • That is, 

I 1 ' 1 

we must have jZ (s)| much less than } c g where 


Z^ ( S ) — hgS + 


Cgs +- 


L^s 

= L s + i 


Li C4B 




14- L^OgS^^ 


a 1 


As « L^Cg , this becomes 

jl-g \ <( 1 


L3 O^uf « 


1 


(5.3.21) 
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L C 
max max 


(L1+L3) (Cg+C^) 




as Cg)) C4 


When ' 


is substituted in h^Q^Lo , we obtain 


K 

3 4 


V4^’kn 


(L^ + L3 )Cg 


«< 1 


as C4«C2 and Lg 4 L^+L^. Thus at 

for *^2 ^4 > i 2 -|_(s) [ <'< “0^g j. Of course for all fO, 

4 

this inequality juay not be true. A.ctually 
corresponds to maximiitti time constant. So the approximated 
system (with 0^ thrown off) response approaches the 
actual system response as the time becomes larger than 
> "the matimum time constant of the system. Note 

' max 

p 

also that L C 1 may be true for some U>> Lu , 

1 'i min 

but not for lJ >>oJ . . So, even for t< 2 < , the 

approximated system response is very close to the exact 
response. When t is very small (just near 0), the 
high frequency modes do contribute to the transient 
response. Thus the approximated system response differs 
much (comparatively) from the exact response, but as 
time increases this difference in the responses decreases. 
Thus we can, as an approximation, neglect O4. when 
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Of course the current source in Figure {5,2,1-c) can be 
transformed into a voltage source to get the network 
for the approximate system as in Figure (5.2,l~b). 

Throughout the discussion in this section, we 
have tacitly avoided the consideration of systems with 
real and complex eigenvalues. However, some discussions 
on such systems will be given in section 5«5. 

5.4 EX/LT\/!PLi! SYSTEK© WITH REAL EIGENVALUES 

Sections 4 and 5 aim at illustrating the 
use of the proposed method. Several example systems 
are considered and for typical cases, responses of 
the actual and approximate systems are plotted in order 
to gain a feel for the extent of approximation made. 

Before proceeding with example systems, let us, 
for convenience, summarize the results stated in section 
5,2. Given a stable system described by 

X = X + f u ( 5 , 4 , 1 ) 

with (A , f) as in (5.2.3). Compute first the elements 
( denoted by b^’ s , i=l ,2, . . . ,n ) of the Schwarz matrix 
( denoted by B ) corresponding to A. Next, carry out 
the comparison process starting with the comparison of 
bn I with bn^g. In general, i-th comparison would 
involve comparing with and this i-th 

comparison is declared successful if the quotient 
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is greater than or equal to k where the 
value of k to be chosen ilepends upon the extent of 
approximation desired* and hence its choice is quite 
subjective. Unless stated otherwise, k = 9.5 is used 
throughout sections 5. 5 and 5,6, A higher value of k 
would imply a higher order for the approximate system 
and hence may result in a better approximation. Whenever 
during the comparison process ^n-i/^n-i-1 iesriS 

than one, terminate the comparison process. Otherwise, 
continue doing the comparisons until ^2 Is compared 
with b^. If i-th comparison is the last successful 
comparison, then infer that the actual system can be 
approximated by a system of order (n-i) . The reduced 
system is described by 

z* = B* z"* + f* u (5.4,2) 

where B* and f* are as given in (5.2,10) through 
(5,2,13), Transform the system in (5.4,2) to one in 
phase-variable canonical form given by 

X* = A* X* + f* u (5,4,3) 

through a transformation (refer, for details, to 
chapter III) , 

z* = P* X* (5.4.4) 

The system in (5,4.3) is the required, simplified system. 

*In a typical dominant-pole approximation method, an 
eigenvalue which is say 25 times larger than the smallest 
eigenvalue (excluding zero eigenvalues) may be neglected. 
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It is wortii mentioning once again that our method 
docs not require tho computdtion ot eigenvalue's "and 
eigenvectors of the actual system. However, to have a 
feel for what is happening in the approximation process, 
in all the example systems, we give first the eigenvalues 
of the original system and then those of the simplified 
system (Note that our approximation or simplification 
method, after eliminating the system ”far-off eigenvalues”, 
modifies the remaining ones to hopefully get a simplified 
system M^ith responses ’’close” to those of the actual 
system) . Comments on the modification (due to our 
simplification method) of the eigenvalues of the system 
are also given for many example systems, finally for 
typical cases, responses of both actual and simplified 
systercsare plotted. 

In this section, we consider a few systems with 
real eigenvalues. The simplification process for the 
first example system is explained elaborately and for 
others, only end results are stated, illso for the first 
example system, the reduction process is performed through 
Davison’s method fssland the results of our method and 
Davison’s method are compared. Towards the end of this 
section, an example system with numerator dynamics is 
discussed. The simplification of this system is carried 
out through three different methods and the results are 
compared. 
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/^^=-0.2, /\2=-0.4, /A2=-0.5, >^=-5.0, >^5=-l0.0 & Ag=“20.0. 

TJie elements bj^’s (1=1 ,2, . . , ,6) of the corresponding 
Schwarz matrix are fotind out as: 

= 0.03590879 

bg = 0.27054737 

bg = 3.18777084 

b = 35.24095232 

4 

b^ = 350.14481408 

5 


36.09999968 
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Now start the comparison process. 

'= 10 
b^/bj = 11 

bg/bg = 11,8 

bgAj^ 1' 7.5 

Note that the first three comparisons are sucoessfiil and 
that the fourth one is unsuccessfiil. Thus, our method 
indicates that the order of the simplified system is 3. 
Using equations (5,2,10) through (5,2.13), the simplified 
system in Schwarz form is written as: 


I 1 

, 

i 2 

; z* 
3 


where 

b* 

1 

'^*3 

and p 



"o 1 o'' 

-b* 0 1 


^1 

"s 


"o ' 

0 

j 

0 -b* -b* 

2 3 

j 

Z* 

i 

^p_ 


0.0359037 

0.27054737 

0.87736430 

0.00078604. 


Apply a transformation 


u(t) (5.4.6) 


-sr* 

^1 


10 0 


rr ^ 

^1 

^2 

; j 

0 1 0 


^2 

^3 ! 


-b* 0 1 

! 

1 lyr^ 

. sj 


■to the system in (5,4.6) to result in a system 


(5,4.7) 
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TJtie responses for \mit step input ^ 2 :^{t), Xg(t}, 
Xg(t) and x|(t), x^(t) and Xg(t) are plotted in 

Figure (5.4,1), Observe that the approxinete responses 
follow closely the actual responses. However, the error 
associated with derivative variables x*{t) and x^{t) 
are considerably large compared to that associated with 
x*{t). This is true for all example systems to be 
discussed, 

Davison’s method LzsJ will now be applied to reduce 
the sixth order system in (5 ,.4,5) to a third order system. 
A computer program has been used to compute the reduced 
system matrices A** and f**. Here, only the final 
results are given: 


A** 


0 10 

0 0 1 

-0.04 -0.58 -1.1 




0.0001137 

-0,0004648 

0,0014694 


and the reduced system is described by: 


X** = A** X** + f** u (5,4,9) 

Note that A** and f** have been calculated from the 
knowledge of eigenvalues and eigenvectors (modal matrix) 
of A, (For details refer to 4** has eigenvalues 

-O.S, -0.4 and -0.5. 
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The approximate system responses x?*(t), x**{t) 
and x**(t) are plotted in the same Figure (5.4,1), 
Corapare the responses of the systems obtained by our 
method and by Davison’s method. x**(t) is closer to 
the actual response x (t) than x^(t;). But x*(t) 
and x*(t) are much better approximations than x**(t) 

O 

and x*'*'(t). It is to be noted that in Davison’s method 
o 

also, error with the derivative variables x^*(t) and 
x^*(t) is considerably large compared to the error with 




(t) . 


Example (5.4.2) ; From now on, only final results will 
be given, Dhstarred quantities refer to actual system 
and the starred ones to approximate system. With this 
notation, all the results will now be presented. 


= -1, /Xg = -2, Ag = -50 and = -50, 

That is, the actual system has a repeated eivenvalue at 
-50 with multiplicity 2. The comparison process resulted 
in 

bg/bg = 38 and 

thus both the comparisons being successful. The reduced 
system has eigenvalues: 


= -1,12617126 

= -1.62795264 

The responses of actual and approximate systems, for a unit 
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step input, are plotted in Pigure {5.4,2). 

Example (5.4,3) ; TJie actual system has 

/^ = -50, = -50, Ag = ~50 and = -100. 

''/e get 

bg/bg =6.7 and bg/b^ = 4.3. 

So neither of the comparisons is successful and hence no 
simplification is possible. This conclusion arrived at 
through our method is consistent with the expected result. 
Example (5.4,4) ; Let us consider a system with 
eigenvalues ’’fairly spread”. 

A^ = -100, Ag = -50, /\g = -30, = -2, /\g = -1 

and /\g = -0.1. 

Vife get in this case 

■^ 5^4 = 9 , b 4 /bg = 33 , bg/bg = 10 and « 80. 

Even though the first comparison is not successful, all 
the other three are successful and hence the simplified 
system is of order two. The eigenvalues of the simplified 
system are approximately -0.084 and -1.6. 

Example (5.4.5) : Let us now consider an example with one 
repeated eigenvalue which is dominant (i.e. to be retained 
in the simplified system). 

A^ = -0.1 , As = -0*1 » ^3 " \ 

A rough calculation indicates that 
bg/bg = 75 and bg/b^ = 65 

i.e. two comparisons are successful. So the simplified 
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system is of order two. 

Let us now consider the last example in this 

section which involves numBrator dynamics. 

Example (5.4.6) f Einally , consider a system with 

numerator dynamics described by 

. (3 4 ) 

TJ(s) (s+1) (s+3) (s+5) (s+10) 


{5*4.10) 


Following the method II given in Appendix II the system 
in (5.4.10) can be equivalently represented by the state 
equations 

h' 

^2 

^3 

^4 

L 

... (5.4.11) 


1 

0 

1 

0 

0 


^1 

1 

0 


0 

0 

1 

0 


^2 i 


0 


0 

0 

0 

1 


X 3 


0 

t 


-150 

-245 

-113 

-19 


X4 


1 


u(t) 


X = 4x^ + Xg 


[4100 


X 3 

x^ 


... (5.4.12) 


where x(t) = L~^[{X(s))j, 

The elements (i=l»S,3,4) for (5,4.11) are 

computed and the comparison process results in 
b g/b g — 8 . S and b g/b — 7.6. 

With k = 9.5, our method shows that no reduction is 


*This example is taken from L38j . The simplification 
is carried out by our method and the response of the 
simplified system is compared with those obtained from 
approximated systems given in [ZS] . 
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possible. However, Meier and Luenberger [_36j have 
approximated the above system by a second order system 
given by: 

X(s) ^ 0.0043 (s~19) (5.4.13) 

U(s) (s+1) (s+3) 

Reference L38j also gives another dominant-pole approximation 
to (5.4,10) as; 

X(s) ^ 0.08 (5.4.14)* 

IJ(s) (s+l)(s+3) 

In order to compare the results of our method with (5.4,13) 
and (5.4,14) , we purposely have chosen a smaller value 
(k = 7,5) for k and carried out the simplification 
process. We get 

X* 0 1 X* 0 

1 1 

= + u{t) 

7* -a| -a* K* p 

and X* = 4x^ ^2 ,,.(5.4.15) 

where 

a* = 3.16302280 

a* = 1.49842202 

and p = 0.00998948 

The eigenvalues of the simplified system are -1.06151139 + 

2 0.57337231, The damping constant of this system is 
approximately 0,89. 

*This is obtained by neglecting s in (s+4), (s+5) 
and (s+10) appearing in (5.4.10), 




TJae responses x(t) (obtained by Ifeier and 
Luenberger’s method), x(t) (a dominant-pole approximation 
of x(t)) and x^(t) (obtained by our method) for step 
input are plotted in Figure (5.4.3). 

It is to be observed that Meier and Luenberger 
assume some knowledge about the input fimction and have 
to solve a set of nonlinear equations to obtain the 
simplified system in (5.4.13). This approach may become 
very difficult if the actual system is of large size, 

5.5 EXA-MPLE SYSTE^^S VflTH REAL AMD COMPLEX EIGENYALEES 

In this section, typical example systems with real 
and complex eigenvalues will be considered and some 
problems associated with simplification of such systems 
will also be discussed. Let us repeat again the following: 
Whereas the proposed method does not require the knowledge 
of the eigenvalues of the actual system, just to gain 
confidence in our method.isnd to gather support for our 
method, we start with system whose eigenvalues are known. 
Having clearly followed the fact as to why we often talk 
about the nature of eJvenvalues, let us understand that 
there are tvi/o distinct cases in complex eigenvalue systems. 
In the first case; the actual system has some complex 
eigenvalues to be retained in the simplified system. In 
the second case, the actual system has some ’’far-off 
complex eigenvalues” v^hich are to be eliminated in arriving 
at the simplified system. In each case, several examples 
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FIGURE (5 4 3) 
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will now be presented. 

Example (5,5,1) ; Consider a system with = -100 and 

~ j 7.75, As expected our method results 

in one successful comparison and the reduced system is 
found to have eigenvalues very close to Xg and As* 

Example (5,5,3) : Consider a system with Ap " -100, 

Ag = -80 and A 35 A 4 = -0.5 ± j 7.75. One may expect 
that reduction by order two may be possible in this case. 

But onr method shows that only the first comparison is 
successfiiL. The reduced system (of order 3) has eigenvalues 
approximately equal to -0.5 + j 7,75 and -44, It may 
be argued that as the real eigenvalue is only about 15 
times (magnitude wise) the complex eigenvalue, further 
reduction is not possible. 

Example (5,5,5) ; Consider a system with 
An = - 1 , A„»A = -0.5 4 - j 3,5, X. = -20 and = -30. 

i- 2 3 ^ 

It is found that 

b^/bjj « 29 , bg/bg tr 6.15 and b^/b^ ^ 1.18. 

Thus, the simplified system has the order equal to 4. The 
simplified system has eigenvalues given by: 

= -1.0311605 , X* = -10.9478410 
and A * A * = -0.5209461 + j 2.5243308. 

The responses of actual and simplified systems are plotted 
in Figure (5,5.1) , 
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In the above example systems, we had dominant 
complex eigenvalues. 

Example (5.5.4): Consider a system with 

« -0.1, Xg = -0.2 and Ag,>.^ * -5 t :J 5. 

The simplified system has order two and it has 

/\!f « -0,10341686 and K* *= -0,18759833, 

4 . 

Thus the eliminated modes correspond to -5 ♦ j 5, as 
expected. Refer to Figure (5.5,2) for responses of the 
actual and the simplified systems, 

- Example (5,5,5) ; Consider a system with 

» -0.5, Xg = - 1,0 and A 3 , A = - 1,0 ± j' 25.0. 

Note that Xg ^ad X^ magnitude wise much larger 
than Ap ^hd Xg* '.though the first compajis on 

is not successful , the second one is and thus the system 
size can be reduced by two. The simplified system has 

* -0,39379313 and = -2.80219074* 

Figure (5,5,3) shows the responses of the abtuai and 
simplified systems, . 

Example {5,5.6) ; Consider a system with 
>1 = -1.6 , , X^ = -0,5 i 3 3.0 aaa = - 0.05 

^ • ' " ' ^ j'O.a 

Hie reduced system has dimensipn 3’ and its eigenvalues are; 
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X* = -1.6834226 and X*> X| = -0.05021498 + 

2 0.19872369. 

Responses of actual and approximate systems are given 
in Figure (5,5 .4) . 

So far in almost all the examples, our method 

uS 

has given expected and reasonable results. Lel^oonslder 
two more examples which present some conceptual difficulties. 
Example (5.5,7) ; Consider 'a system with 

^2 — —1.0, s yKg ~ —0,5 + 3 3.0 and *=-8,0 + J 0,2, 

Note that A 4 and are, magnitude wise, not much larger 

5 

than/Xg and However, our method results in 

b^/bg = 11, ^^ 3/^2 “ ^ad i>2/^i " 1.13. 

Thus our method says that the given system can be simplified 
and the simplified system has order 4. The eigenvalues of 
the simplified system are found to be: 

^\* = -1.1467420 , = -2.9441188 

and 

= -0.44205414 + 3 3.09066888. 

Note that A* is close to A 2 -^3 ^^'3, A * are 

close to An ^^'3 Ar,* the following comments are 

now in order: 

1 . Sven though the 'magnitude criterion' (which is found 
to be successful in all the previous examples) sx:iggests 
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.that no reduction is possible in this case, o\xr 
method does indicate the possibility of reducing 
the system by one order. 

2, In the simplified system, one eigenvalue in the 

complex pair and /\g is eliminated and the 

other modified drastically. 

However, the responses shown in Higure (5,5,5) 

show that the approximation made is fairly good. It 

may be argued that the effect of the eigenvalues 

and (This pair has a damping constant very close 

5 

to 1 and can be thus regarded as approximately equivalent 
to a real repeated eigenvalue at -8,0) is approximately 
realized by having a single eigenvalue at -S. 9441 188, 

A rigorous justification of this statement needs further 
investigation. 

It is interesting to ask as to what will happen if 
the damping constant (f associated with and Ag is 

decreased from its present value. In all cases with 
6 = 0.85, 0.80, 0.707 and 0,6, it is found that only 
one comparison is successful thus resxlLting in the same 
situation. When 6 = 0.15 (the imaginary part of 

and Ag = 50), our method indicates that the size 
of the system can be reduced by two. Observe that in 
this case, the magnitude of "^5 much larger 

than those of As* 
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EXACT AND APPROXMATE RESPONSE FOR EXAMPLE (S 5 V) 


FIGURE (5 5 5) 
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.Example (5»5,8) - Consider a system witJi 

/'!’ Az ~ ±53 and Ag, X = -8,0 + J 0,2. 

This also resulted in one successful comparison and the 
simplified system has 

= -3.6951548 and = -0.4117478 + J 3.04310048. 

A-11 the comments given in Example (5,5.7) may be made for 
this Example also. Responses of actual and simplified 
systems are given in Eigure (5.5.6). 

Before ending this section, let us observe the 
following: From an the examples presented in section 5.4 
and 5,5, it is inferred that given any system described 
as in (5.2.3), without v/orrying about the eigenvalues, 
the system may be simplified by the proposed method and 
that the simplified system response is ’’close” to the 
actual system response. However, the error associated 
with derivative variables (x^, Xg, etc.) is considerably 
large compared to that associated with the variable x*. 

Results of the typical examples considered in 
sections 5,4 and 5.5 are summarized in Table (5.5,1). 
The various responses for actual and approximate systems 
for these examples were computed using the program 
’’Modified Digital Analog Simulator” written for IBM 7044, 
For simulation of the system in Schwarz form, refer to 
section 2.3. 
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TABLE (5.5,1) 

ACTmi AM) APPRO XimTE SYSTEMS DATA 


Example 

Actual System Data 


Simplified System Data 

Fo. 

Or- 

der 

Schwarz Eigen- 

elements values 

Or- 

der 

Schwarz 

elements 

Eigen- 

values 

(5.4.1) 

6 

0.03590879 -0.2 

3 

0.03590879 

-0.16883593 



0.27054737 -0.4 


0.27054737 

--0,35450246 + 
50.24782361 



3.18777084 -0.5 


0.87736430 




35.24095232 -5.0 






350.14481408 -10.0 






36.09999968 -20.0 




(5.4.2) 

4 

1.83335350 -1.0 

2 

1.83335350 

• 

-1.12617126 



72.92392704 -2.0 


2.75412392 

-1.62795264 



2727.24269056 -50.0 






103.00000000 -50.0 




(5.4.6) 

4 

1.50000000 -1.0 

2 

1.50000000 

-1.08151139 + 
jO. 57337231 



11.40000000 -3.0 


2.16302280 




100.10000000 -5.0 






19.00000000 -10,0 




(5.5.5) 

5 

3.20538708 -1.0 

4 

3.20538708 

-1,03116130 



3.80996440 -O.S+j 

2,5 

3.80996440 

-10.94784100 



23.39810976 -20.0 


23.39810976 

-0.52094610 + 
52.52433080 



677.08653056 -30.0 


13,02089472 




52.00000000 
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119 


(5.5.4) 

4 

0.01940083 

-0.1 2 

0.01940083 

-0.10341686 




1.45632678 

-0.2 

0.29101519 

-0.18759833 




51.54427072 

-5.0+j 5,0 






10.29999992 





(5.5.5) 

4 

0.86720758 

-0,5 2 

0.86720758 

-0.58379513 




267.70421504 

-1.0 

2.59598388 

-2.20219074 




360.92857344 

-1.0+j25.0 






3.50000000 





(5.5.6) 

5 

0.03964259 

-1.0 3 

0.03964259 

-1.68342260 




0.17143613 

4.72225436 

5.55916672 

2.09999996 

-0.5+j 3.0 

-0.05+00.2 

0.17143613 

1.78385260 

-0.05021498 
jO. 19872369 


(5.5.7) 

5 

4.05029576 

-1,0 4 

4.05029576 

-1.14674200 




4.56506204 

-0.5to3.0 

4.56506204 

-2.94411880 




8.12519726 

-8.0+j0;2 

8.12519736 

-0144205414 
j 3. 09 06 6888 




89.54944384 


4.97499080 





18.00000000 





(5.5.8) 

4 

7.71143768 

-0.5+j 3.0 3 

7,71143768 

-3.69515480 




4.76150336 

-8.0+0 0,2 

4.76150336 

-0.41174780 
j 3. 043 10048 

t 



76.81705856 

17.00000000 


4.51865044 




]v}OTB : 

TJie Schwarz elements stated in the above Table are in the 
order bi,bp...,b for actual system and bj,b*,...,b*^ for the 
simplified system inhere n and a* represent Respectively the 
order of the actual system and that of the simplified system. 
Except for Example (5.4.6), k = 9,5 was used (See section 5.4) 
for all Examples and for Example (5.4.6), k = 7.5 waia used. 
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5»6 PROBL®© FOR FOTRRl INVESTI&A^TION 

A few problems associated with the, proposed 

M. 

simplification method may be investigated further. The first 
one is a systematic extension, of this method to overcome 
some of the difficulties encountered with complex 
eigenvalue systems. Second, the method may be suitably 
modified to bring multivariable systems also within its 
range of application. In this connection, the Schwarz 
canonical form for miiLtivariable systems and its 
network interpretation discussed in chapter IV may be 
useful. Third, procedures may be developed to estimate 
the amount of error between the exact and approximate 
responses. Finally, it may be possible to determine 
a performance index which becomes optimum for the proposed 
simplification method. 

5.7 CONCLD3ION 

The proposed method has the major advantage that the 
system eigenvalues and eigenvectors are not to be evaluated. 
Schwarz canonical form has been used both for determining 
whether the system can be simplified or not and for 
simplifying the system. Schwarz form is also useful in 
detecting redundant state variables and eliminating them, 
(see Appendix V), 

The simplified system has a few spectral energy 
density moments exactly same as those of the actual systems. 
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For details, refer to [42, 25] . 

Several oxamples iiave been provided, to illustrate 
the simplification method. However, it is to be emphasized 
that as is common with many simplification procedures, 
the proposed method also does not give reasonable results 
for all inputs. 
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CHAPTI3R - YI 

^ CAnOlJICiiL FORM FOR LINEAR Tllffi- INVARIANT 
DI3CRETE-TII\CE SYSTEM®! 

6.1 INIRODDGTION 

Th-: Schvfar^ canoT-ical form for linessy time -invariant 
continuous systems Las mainly the following three properties: 
(1) stability of the system is inferred by inspection of the 
elements of the canonical matrix; (2) Liapunov function for 
the system is constructed quite simply snd (3) the elements 
of the canonical matrix can be computed from the elements 
of the first column of the Routh array . It is the aim of 
this chapter to develop- quite analogous to the development 
in the continuous case- a canonical form for linear time- 
invariant discrete-time systems, in fact, the canonical 
form proposed in this chapter has the first two properties 
stated above. Analogous to the aforementioned property (3), 
attempt is made to compute the elements of the proposed 
canonical matrix from those of the Jury table. Presently 
this attempt is not completely successful. However, it is 
hoped that in future such a linh- between the proposed 
canonical matrix and the Jury table- would be established. 

Contents in this chapter will now be rapidly stated. 
Section 2 sj'stemat j cally derives the canonical form. It 
also discusses some properties of the canonical matrix. 

Use of Ghis canonical form in investigating the system* 

''‘For couvonionce, instead of every time saying 'linear 
time-invariant discrete-tirao system*, the word 'system is 
used tb-r-oughoub thi s chapter. 
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stability as well ag in constructing a Liapunov function 
for the system is discussed in section 3. Section 4 gives 
the relationships between the elements of the canonical 
matrix. The fairly lengthy section 5 enlists a few 
methods of computing the elements of the canonical matrix. 
Finally, a few problems are posed for further investigation 
in section 6. 

6.2 PROPOSED CiiiroiTECAL FORM 

The treatment given in this section closely parallels 

1 r- 

'_45j . The concept of matrix bilinear transformation 1.45, 
47 J is freely used in the sequel. To make matters simple 
in the actual derivation of the canonical form, the 
properties of the matrix bilinear transformation will first 
be enunciated without proofs (For proofs, refer to (_45J ). 

V is defined ag the matrix bilinear transformation 
of U if 

Y = (p'+l)(U-T)“^ (6.2.1) 

ITote carefully that the above definition presumes that 
■ 5)"^ exists. In other words, the matrix p" does 
not have an eigenvalue (or multiple eigenvalues) at 1*. 

The statement of the properties of the matrix bilinear 
transformation is now in order, 

^Throughout the discussion in this chapter, it is 
assumed that the discrete system under investigation has 
no eigenvalues at 1. For a stable system, this is 
automatically satisfied, as a stable system, if continuous, 
has all eigenvalues in the left half -plane and if discrete, 
within the unit circle. 




1 


coimaute. Thus 


. (U + I) and (U - 1)“^ 

I = (u + I) (y - i)~i = (y - 1)“^ (y + i) (6.2.«) 

2. A raatrlx commutes with its bilinear transform. That is, 

yv = TU (6.2,3) 

3. The bilinear transformation possesses the property 
of reciprocity. This means that if V the 
bilinear transform of y , then y is the bilinear 
transform of v. That is, given (6.2.1), 

y - { Y + IMY - I)~^ (6.2.4) 

4. If 

U = S T , (6.2,5) 

then (6,2.1) implies that 

Y = 3 (T I) (T - l)-l S“^ (6.2.6) 

5. If all eigenvalues of a matrix are in the left half 
plane (unit circle), those of its bilinear transform 
are in the unit circle (left half plane), 

6. A matrix and its bilinear transform ha^the same 
eigenvectors. Further, the bilinear transform of 
a non-derogatory (derogatory) matrix is itself 
non-derogatory (derogatory)’^. 

*A matrix is said to be non- derogatory, if it has only 
one eigenvector per eigenvalue and derogatory, otherwise. If 
a matrix is non-derogatory, it is similar to its companion 
form (phase-variable form). 
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The equation (6.2.1) can be put in a more useful 
form as under- 

2(u - l)-l (6^2.7) 

Before proceeding to the derivation of the canonical 
form, a clear statement as to v/hat is aimed at may prove to 
be useftiL. Given any stable system (discrete) described by 

1 ~ ^ — it ^ 6 # 2 , 8 ) 

where A is non-derogatory . ’'Ve wish to transform (6,2,8) 
to a canonical form 

where W is similar to a. m what follows, we consider 
W as the bilinear transform of the Schwarz matrix B of 
the continuous system obtained from (6.2,8) through the 
bilinear transformation. In the next section, we shall 
see vjhat the stability of (6.2,8) or (6.2.9) implies in 
terms of the elements of W. 

Define 0^ as the companion form of any matrix 
(non-derogatory) S. We are now ready for deriving the 
canonical form. As a is non-derogatory, 

= N A (6.2.10) 

Let D be bhe bilinear transform of 0^. Then by the 
property 4 stated above, 

D = ^ (A + I) (A - I)"^ 


( 6 . 2 . 11 ) 
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By property 6, it is clear tiaat D is also non-derogatory. 
Therefore , 

5d " 

i.e. 

—1 1 1 

Sd = M N {A + I) (A - I) (6,2.12) 

Transform to its Schwarz form B. That is, 

B = P (6.8.13) 

or 

B = (P M N) (A + I) (A - I)-^ (P M N)“^ 

Denote 

S * ? M N (6,2.14) 

Then 

B = S {^1 + 1) (^_^ - I)"^ S"^ 

Invoking the properties 5 and 4 , 

(B +1) (B - 1)“^ = S A S“^ 


Define 

W = (B + I) (B - I)"^ 

Thus 

W = S S"^ 


( 6 . 2 . 15 ) 


(6.2,16) 


Thus we have shown that W as defined in (6,2.15) is 

■li. 


similar to 
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Mow a few comments are in order. Let G- be the 
bilinear transform of A. Thus, 

G = + I) (A - I)"^ (6.2.17) 

Transform G to its Schwarz form B. 

§ = ^ 1 (A +1) (A - 1)“^ H (6.2.18) 

Therefore, from (6.2.15) 

W = A H (6.2.19) 

Mote that in equations (6.2,17) through (6.2.19), 

transformation to companion form is not involved. One may 

wonder whether non-derogatorihess restriction can be relaxed 

or not. It cannot be relaxed because exists only 

when G (hence A, by property 6) is non-derogatory. 

However, a problem where a is derogatory can be tackled. 

for details refer to [ 4 : 5 \ . 

The several steps in the derivation of the canonical 

matrix W may be well illustrated with the help of the 

Figure (6.2.1). The scalar bilinear transformation 
s+1 

(z = — r) which is used to find the elements of from 

those of Cg^ becomes very tedious when the order of the 
system becomes large. However, an improved technique 
which has been recently proposed by power [47j may be used, 
to achieve the same purpos'h, A method of computing D 
from C is given in section 6.5'. The transformation 

— Q 



SCALAR bilinear TRANSFORMATION 



transformation to the proposed canonical, form 

FIGURE ('6 2 1; 
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matrix p wliicli transforms to B is discussed in 

ciiapter III. The elements of W may he computed 

from the elements h^^’s of B in two ways. Method 1 uses 

a set of rules which help in computing w, .’s from b^’s 

ij 

and this is discussed in section 6,3. Method 2, which makes 
use of recursive relations between w, -*s and bi ’s is 
presented in section 6^5, In section 6,5, attempts are 
made to express 'bj_'s directly in terms of the elements 
of (say aj_’s} and to relate with through 

the elements of the Jury table of the characteristic 
polynomial associated with C^. 

Some interesting properties of the canonical matrix 


W 

u 


(for a stable system) are stated below: 



•» g '-.J 


( 6 . 2 . 20 ) 



i> j 


or ■ 

L k-j 



( 6 . 2 . 21 ) 



( 6 . 2 . 22 ) 


(6.2.23) 
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Properties S and 3 will now be explained in words. All 
the superdiagonal elejnents of the canonical matrix are 
negative. Elements of the canonical matrix lying in odd- 
ordered subdiagonajs are positive and those lying in even- 
ordered subdiagonals are negative. While dealing with 
many numerical examples, it is found that 

<.^1 » Wj )> 4 1 {6.2.24) 

where and w^ are respectively the i-th and the 

j-th columns of W. Proof of (6.2.84) in general is 
not available presently. However, in what follows proofs 
of (6.2.20), ( 6 , 2 , 32 ) and (6.2,23) are given. 

Proof of (6.2.20) : 

The commutative property given in equation (6.2.3) 
of the matrix bilinear transformation is used below in 
proving (6.2,80) by induction. From (6.2.3), 

WB = BW (6.2.85) 

Equating the corresponding elements of W B and B W 
on the leading diagonal gives relations between elements 
of W on the first sub- and superdiagonals; equating 
corresponding elements of W B and B W on the first 
sub- and superdiagonals gives relations between elements 
of W on the second sub- and superdiagonals; and so on, 

( W B ) . ^ = ( B W ) . . 

- - 
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gives 




' 

• • ♦ ( 6 3 *26 ) 

(6.2,26} implies 

- ->>3 »j(34.i), 

(6.2.27) 

provided 

W -b W 

d(j-l} M (j-Dj 

(6.2,28) 

Now, (W B}^^ = (B results in 

Wgi » Thus 


from (6.2,27), Wgg = -bg Wg 3 and thus W 43 = -bg W 32 
Q'lid so on. So by induction it has been shown that 


;j+'l) j‘- ~ ''*'^^^5+1) * {5=i»2,.,,,n~l) 

.»• (6.2,89} 




Next, from 

we get 

"^0(1-2) -^"(j+DCM) “ “^j-i'"(j-.i)(j-i) 

... (6.2.30) 

Similarly 


(W B),, . . = (B W) 

- - ( J-l)l - - (3-1) J 


yields 




{6.8.31} 



13E 


Miatiplying (S.S^Sl) by b 

‘*’3-1 “'(a-DU-i) 


and rearranging' result in 

"^( 3 - 2)3 

+ b . , w . . 

3-1 33 


Note tiiat the right hand sides of (6.2.30) and 
are identical and therefore ■we have: 


^(3+l)(3-l) ~ ^3 ^3-1 "'(J-Dii+l) 


(6.2.33) 


provided 

^3(3-2) = ’’^3-1 ^3-2 ^(3-2)3’^^^ (6.2.34) 

But from (W B)g^ = ^-?^12 " 

Qfter some manipulation, •we get 


W31 = 1>i^2^13* 

Therefore from (6.2.33) and (6.2.34), 'w^g = hgbgWg^ 
which in tiirn Implies w^^ = b^b^w^^^ and so on. So, 

by induction it has been shown that 


w 


(3+l)(3-l) “ 

... (6.2.36) 


The same procedure may be continued for the second sub- 
and superdiagonal elements, third sub- and superdiagonal 
elements and so on to obtain the ecuations 
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w 


( s+i) s 


(-b ) (-b )...(-.b . .) w . 

1 1+1 s+i-l i{s+i) 


1=1, S, . . . ,n-l 

s=l ,2 , . . . ,n-i (6,2.36) 


(6.2,36) can be rewritten as in (6,2.80). Tills completes 
the proof. 

Proof of (6,2,28) - 

Define F = (B - Then 


F (B - I) = I 


(6.2.37) 


(6.2.37) implies [F(B - I))^^ = 0 , 1 < n. That is 


f . , , . - (1 + b ) f . = 0 

i(n--l) n in 


or 


h(n-l) = (1 + V ^ 


in 


(6.2.38) 


Equating I F(B - I) I . , \ to zero, 

- - - - ivn-1 > 

^l(n-2) ” ^i(n-l)‘^ "^n-l ^in 


v;e get 


(6.2.39) 


Repeating this procedure, we get 
^i(n-k) “ ^i(n-k+l)'*’ "^(n-k+l) ^i(n-k+2)’ 


(k=2,3, ... ,n-i-l) 
... (6.2.40) 


From the above equations (6.2,38) - (6.2,40), it is clear 
that all the superdiagonal elements lying in the i-th row 
of F have the same sign as that of fj_^. As W is the 
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matrix bilinear transform of B, from (6.2.7) we get 

W = I + 2(B - I)-l (6.2.41) 

Therefore (j > i) = 2 f^^. Thus proving (6.2,22) 

reduces to showing that 0 (J ^ i) hence it has 

to be shown that f 0. 

Referring to Appendix I Eq.ua tions (A. 1.3) , 

= (1/B(n)) F(j-l) (-1)^'“^^ (6.2.48) 

AS Bin) has the sign of (-1)^ and F(i-1) that of 
(-l)j-l^ it is readily seen from (6.8.42) that f C 0 » 


This property of W given in (6,2,23) is obtained 
by combining (6.2.21) and (6,2,22). From (6.2.21), 


(s+i)s 


s(s+i) 


(s+i-1) 

7T 

k=s 


(6,2 ,43 ) 


From (6.2,22), Therefore, for a stable 

system (bj^’s>0), '’the’sign of is that of (-1)®“^!. 

Hence the relations (6,2.83). 


6.3 STABILITF IWESTIGATION 

In this section, first, simple rifles to compute 
the elements of W from those of B will be stated . 

Second, the proposed canonical form is used in constructing 
Liapunov functions for the system, in section 6.2, a stable 
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Rule I To eval*uat8 £\ 

— jQ 




n 


1+ 


n 


Sbi 

1=1 


n-2 n 

^21 s: Vj 

1=1 j=l+2 


n-4 n-2 n 


* X ZI ZI 

i=l o=i+2 k=j+2 


•f 


... (6.3.2) 


For example, •with, n = 5, the abo’ve rule gives 


^ - 1 + h + b + b + h + h + b b_ + b.b^ + b,b^ 

^ 131415 

+ bgb^ + bgbg + bgbg ^ T^i^gb^. 

We observe that can also be found recursively by 

using the relations (Appendix I). ^ 

B(0) = 1 

B{1) = -(1 + h^) 

B(k) = -B(k-l) 4 B(k-2), k = 2,5,...,n^, 


> 


(6.3.3) 


Then 

A = B(n) 
n 


However, we require to be written explicitly as in 

(6.3.2) for applying the other rules to be stated. 

Rule II : To determine the diagonal elements of W 
From (6.3.1) , 


w ’ . 
ii 


w 


11 


A 


n 


Replace b^ by (-bj_) in the expression for 




to get Replace tg by (-bg) in w’^ to get w^g. 

To determine w' (k = 3,4,...,n), replace b ^_2 and b^^ 
by (-bj._ 2 ) and (-bj^) respectively (k-1 ) * 


Rule III To determine tbe superdiagonal elements of W 
^lu+l) is obtained by replacing 

^i+1’ ^i 1 (when i=l, these are bg and b^ 

only) by zeros in the expression for and multiplying 

the result by two. Then, any wl . , jJ>(i+l)j i=l,2,...,n, 

X J 

j = i+2, (i+3) , . . . ,n is found out by Just replacing b^ 
by zero in the expression for w* 

i(j-l) 

Rule lY; To determine the subdiagonal elements pf_ f 


i-1 

W (i> J) = /X (-V 

■' k-3 


A few comments are now in order# First, for 

applying the rules II and III, we must explicitly write 

A- in terms of bi’s as shown in (6.3.2). Second, 
n 1 

these rules are quite simple both for remembering and 
for applying. Finally, it is possible to derive relations 
which help us to algorithmically determine w^^’s from 
the knowledge of bj^'s. This will be done in section 6.5. 
Let us now illustrate the application of the above rules 
by means of a simple example. 






Given 


= 1, bg = 2, bg = 3 and b^ = 4, 

I'D compute W. 

Now, applying the above rules, we obtain 


A 

“"4 

= 

1 + + tg + bg 

+ b. + b-b_ + b^b. + b^b. = 26 
4 13 14 2 4 

w* 

11 

= 

-1 + bj -bg -bg ■ 

-b. + bib„ + bib^ -bnb. =-10 

4 i O 14 <0 4 

W’ 

22 

= 

““1 *^*^2 3 

-b^ + bibs + bib4 + b2b4 = 10 

"^33 

=r 

-1 -bj * bg . bg 

-b^ + b^bg -b^b^ + bgb^ = 6 

^44 

= 

-1 -bj -bg . bj 

+ b^ + b,b„ + b-b. + b„b. = 18 
4 13 14 2 4 

w ’ 

12 

s= 

S( -1 -bg -b^ ) 

= -16 

wl„ 

23 

= 

8( -1 -b^ ) 

= -10 

w* 

34 

= 

S( -1 -bj ) 

= -4 

w’ 

13 

=: 

2( -1 -b^ ) 

= -10 

^14 

t- 

2( -1 ) 

= -2 


= 

2( -1 ) 

= -2 


= 

-'^1 “Ie 

= 16 



'’sh ™l3 

= -20 


= 

■'’s ’'ss 

= 20 

ki 

=r 

-h'=s'’3 "h 

= 12 

^42 


'"E4 

= -12 

^'45 


-’’3 ’’k 

= 12 


12 



inus 



-10 

-16 

-10 

-2 

16 

10 

-10 

-2 

-20 

20 

6 

-4 

IS 

-12 

12 

18 


Observe that all the superdiagonal elements of W are 
negative. That all elements on the first and the third 
subdiagonals are positive and that those on the second 
subdiagonal are negative. Also 


<^!!i . H2> 

■ 'Szy 


334 

36z36 


< 1 


36 

36x36 




1 


’ " 4 > 

<*S ’ !!4> 
’ ™4> 


384 

36x36 

-5-6 

26x26 

284 

26x26 


1 

1 

1 


253 < 

26x26 


Thus we see that ("w. ,w\<l,i7^3. In fact, 
when i = j = 1, this inner product becomes 900/676 > 1. 
Let us next consider the problem of constructing Liapunov 
functions for the system. 

Consider the discrete system described by {6.2.8), 



140 


Transform tiie system through _ - j 

y, = S X, (6.3.4) 

-k - -k 

so that the equation (6.2^9) represents the transformed 
system, Le’|; the Liapunov function for the canonical 
system i^ (6^2^9) he giv^nh^: 

\ “ (Xic ■ i !k> 

where the form of L will he chosen later. The Liapunov 
function expressed in terms of the original system 

state vector becomes 

The first difference of the Liapunov function may he 
readily written as 

Vi ' (h ’ 

or 

\+l ■ \ ' . s'' L W - L) 3 ;,£> (6.S.8) 

We note that L in (6.3.5) has to be appropriately chosen 

such that the first difference of the Liapunov function 

given in (6.3.7) becomes negative definite (or negative 

semidefinite with the additional qualification that 

Y -T does not vanish along any solution trajectory 
k+i k 

of (6,2,9) except the origin). Choose tentatively, 



i^fci 


aiag. (b^b^.-.b^l b^_^b^, b^l 

... (6.3.9) 


Denote 


\r L w 


'j\r - L = ^ 


(6.3.10) 


Then we can write that 


Vi - \ 


^kl " ^kS ■*■ ^”^3^ 


-(-%)* yi,4 - 


... (6.3.11) 


As will he seen, , qgg etc. are negative and thus 

the square roots of etc, are all real 

nuflihers. The expression in (6.3,11) is negative semidef inite 


It may he verified that 




i - 1,2, ... ,n 

h = 0 if k 0 
k 


(6.3.12) 




i? i 


(i+j+1) 




(i = l,2,...,n-l) 

(3 = i+1 , i+2,..,,n) 


(6.3.13) 


q = q , (j = 1 ,2 , . . . , n-1 ) , (i - J+l> 0 + 2 ,,. ,,n) 
iO O'i 

. . . ... (6.3.14) 

X / J 



It can be shown that 


14:3 


(6,3,15) 

Owing to (6.3,15), only the first principal minor of Q 

is different from zero. This principal minor, is 

negative in view of (6,3.12). Therefore, Q is negative 

sem.ld.ef inite . To illustrate whatever has been said so 

far, let us consider a simple example. 

Exam£^; For n = 3, construct ^ and Q, y . 

-k ~ -k 


Applying the formulas given in (6.3,12) - (6.3.14:), 


-(Sb^bgbg) (Sb^bgbg) 

-(Sbjbgbj) 


(Sbgbg) 

iBbgd+bj^)! 

il 

1 

(Sbibgbj) -(Sbgbg)® 

(2bgb3) 

4 

(®V3> i 

^Sbjdrtj)] 

f 1 ^ 


-(Sb b b ) (2b„b„) 

isbjd.bijJ- 

where A 

0 

= ^ ■*’ ^2 ^3 ^1^3 



Uote that ^2 2 ~ ”^1 ^32 and *13 p ~ ”""^1*^32* Now, 

from (6.3.11) and (6.3,12) 


V 


k+1 


? 


k 


T 

y ^ y 

-k - -k 




j 2b2_b2h3 


■i . 


^ A 3 


kl 


2b2b3 

A 3 


yk2 


2b 3 ( 1+b ) 

A 3 


<L 


k3 



Ib is yet to be concluded, that the negative 
seniaef initeness of Q assures the asymptotic stability 
of the system in (6.2,9). Following reference [4:6"] , it 
may be po-sible to give a general proof for the fact that 
the first difference ® identically zero 

for any general sequence of vectors. Hoivever, in this 
section, a simple example will be considered to illustrate 
this idea. 

Example : Consider a system with n = 3. Let b^ = 2, 

"'^2 ~ ^ '^3 ~ system in the canonical form 

can be written as 


^(k+l)l 


2/18 -10/18 -2/18 


^kl 

^(k+l)2 


20/18 8/18 -2/18 


^k2 

./(k+l)3_ 


-12/18 6/18 12/18 




and, 

1 2 

Tr Tr ) ^8 24 24 L 

- Vk = -{« hi - li haj 

The veetc 6' (see reference ) is given by 

T ^ ( 48 , _ 24 , 24 ) 

18 18 18 

As in [46j , let us form a matrix D as 



That is 


i44 


D 


24/9 

-56/87 

-568/243 


- 12/9 v 

- 44/27 
-140/243 


12/9 . 
12/27 
172/243 


We find that |d| / o. So following [46] , we oonolnde 

that identically zero for any general 

sequence of vectors. 

It may he of interest to express L and Q in 
terms of themselves. From the relation (6.2,20), 

the following may be readily obtained; 

m 

b^ = (-1)^ , m£(n-l) (6.3.16) 


k=l ^l(m+l) 

The elements of the diagonal matrix ^ are given by 

-^ii “ ^n^n-l‘**^l 


or 


i 


li 


Vn-r*^^8^1 

b b b . , ,b 
12 3 i-1 


(6.3.17) 


From (6,3,17) and (6,3.16) we get 




ii 


= (- 1 ) 


n-1 


!!nl 


w 


li 


n 


w w 
In il 


(6.3.18) 


Using (6.3.12), (6.3,16) and (6.2.20), we can obtain 


kk 



"^nl^lk 

j. ’"81 


"^41 

1 

11 

2(-l)^-'i b^ ^ 

^ w- w, ■ 
L In kl 

^ ' ’"is 

' ’"is "si 

w 

14 


"^(k-Dl ^l(k-8) 


w 


l(k-l) ^(k-2)l 


. . . (6.3.19) 
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Froia the definition of W , it may he easily seen that 

V/ B = W + B + I (6.3.80) 


Equating the nn element on both sides, we get after 
rearrangement 



Thus equations (6.3.18) - (6.3.21) with (6.3.13), (6.3.14) 

and (6.3.11) determine L , Q and y'^ Q, y in terms of the 

- - -k - -k 

elements of W, 

Referring to (6.3.6), we note that to construct 
the Liapunov function for the system in (6.2.8), knowledge 
of S , the transformation matrix is required. From (6>.2.14), 


S = P M N 

The determination of M and iJ are discussed in * 

For determining M , D given in (6.2,11) has to be first 
evaluated. The transformation matrix P has been 
discussed in chapter III. Reference [2oJ has given a 
method to compute a number of transformation matrices M 
and N. Hence a number of Liapunov functions can 

be constructed. 

Suppose W has been computed somehow. In this 
case, instead of finding P , M and N in order to 
evaluate S , another method is possible. If 



146 


2a = N A If^ 


transformation T can be found such that 


Cg = 2 W T"^ 


Then 


S = T~1 N 


(6.S.S2) 


Finally, let us derive a set of conditions on the 
elements of W for the asymptotic stability of the system 
(6.2.8) or (6.2.9). It is known (chapter II) that the 
system (discrete) is asymptotically stable if all 
(i=l,2, , . , ,n) . From (6.3.16) 


w 


nl 


w 


In 




n 


) (- 1 )^“^ ( 6 . 3 . 23 ) 


That is 




(6.5.24) 


Now ^ 0 implies that 

(-1)^ * b / 0 

nwi^ Wii 

From (6.3.24) this in turn implies that 


(-l)iV (Wji/w^^) > 0 


As ^^'- 11 ^ want this to be so), the above gives 
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(-1) > 0 , i=2,3,...,ii (6.3.25) 

Frora (6.3.84) 

'^11 " (-1)^**"^ (6.3.26) 

As we have forced ["by (6.3.25)J (-l)^'*’! (w, /w . ) > 0 , 

in (6.3.86) ^ 0 if That is, [see equation 

(6. 3. El)] 

nln-i) 

:; > 0 (6.3.27} 

1 - 

m 

EQuations (6.3.85) and (6.3,27) represent the conditions 
on the elements of W for the asymptotic stability of 
the system. 


6.4 RELATTOIB BETWEEN THE ELEMENTS OF W 

From the restilts of the previous section, it may be 
observed i-hat the n^ elements of W are uniquely 
determined in terms of the n elements b^ (i=l ,2 , , . . ,n) . 

O 

So it is reasonable to expect that only n of the n 
elements (i,j = l,2,..,,n) are independent and that 

the other elements are related to these ’’independent” 
elements. Already we have seen that the subdiagonal 
elements are related to the superdiagonal ones through 
bjj^’s as in (6.2.20). In this section, we are going to 
state the relations between subdiagonal elements and 
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between superdiagonal elements of W. Relations between 
tile elements on the tri-diagonals of W (leading, first 
sub- and first super-diagonals) will also be given. linally 
it is illustrated that the knowledge of the first row 
elements of W only will suffice for calculating the rest 
of the elements of W. A simple example is also provided 
at the end of this section. 

The main equations used for arriving at the several 
expressions to be given later are: 

Vf B = W + B + I (6.4.1) 

B W = B+ W + I (6.4.2) 


The derivation of the results of this section involves 
straightforward (though tedious) procedure of equating 
the corresponding elements on both sides of (6.4.1) as 
well as those of (6.4.2). Hence, we will state the results 
only. Of course, the equations used for obtaining such 
results will be mentioned. 

First, we shall write down the relations obtained 
by equating the first row and the n-th row elements of 


(B ^^) and (\^ ^ B + p . 


•». 



V/ 

21 

1 

w 

11 

■Wor, 

2S 

1 — 1 

w 

12 

^2k = 

w 

Ik 

, k = 3,4,.,.,n 


(6.4.3) 
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"n-i 


1 . -u \ 

“ “(u-Dk ’ 


nla-l 

1 - w. 


1 ^ 

n-1) (n-l) 


:6.4.4} 


^ ■*■ ^rtn **■ ^y, ("W - 1) 

b = b nn ' 

a-1 ~~~ T 

(n-l)n y 

E'luatiiv, tlie elements on tJie tridiagonals of both sides 
of (6.4.S) resiiLts in 


[k+l)k ~ ^(k+8)k 
1 - w. 


, (k = l,3,...,n-2) (6.4.5) 


k(k+l) 


, (k=l ,2, ... ,n-2) 


... (6,4,6) 


c+a k+J 


c+1 k+J 


k(k+2) 


■, (k=l,2, . . . ,n-2) 


... (6.4.7) 

Hote that equations (6.4.5) - (6.4.7) express tlie relations 
between the elements on the tridiagonals of W, We shall 
now state the relations between the superdiagonal elements 


of 


^(k-f3)m ~ ^(k+l)m 


(6.4.8) 


( k = 1,2,..., n-3 ) 

(m = k+3, k+4, . . . ,n) 
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(6.^.8) is obtained by equating the elements lying above 
the first superdiagonal (except those in the first row) 
of (3 W) ana (b + W + I) . Relations (6.4.8) can be further 
manipulated to obtain the following simple and elegant 
relati ons; 


... .ha (e.4, 

’"ki+D ’* 1(1+2) ’"in 

i = n“2 

Let us explain what (6.4.9) means. The superdiagonal 
elements of ¥ have the interesting property that each 
element in any i-th row (i ^fn-1) bears the same ratio 
with the element in its own column in the first row as 
any other superdiagonal element in the same i-th row 
does v/ith its corresponding first row element. Next 
lot us demonstrate that a similar relation exists among 
the subdiagonal elements of W. 


of 


Equating the elements below the first subdiagonal 
(V-; B) and (B + W + I) [see equation (6.4. 1)J yields 
w 

^ I5; = 3,4, ... ,n 

^k2 



^m(k-l) - ^mk 
^m(k+l ) 


s. 


k = 2,3,...,n-S 


/ 


m = k+2, k+5,...,n 


(6,4.10) 
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I'li.'-tioiio can be further manipulated to obtain 

the simpit' relations 


1 + 1)1 


"^(1+2)1 ^ ^ ^ni 

’^( 1 + 2)1 %1 
^ = 2,3, . ,j2_2 


(6.4.11) 


Hlqur^tion (G.4.11) means that the subdiagonal elements of 
W ha 10 the interesting property that eich element in any 
i-th o^lunin (j <n-l) bears the same ratio with the element 
in itr. own row in the first coliimn as any other subdiagonal 
element in the same j-th column does with its corresponding 
first column element. 

Suppose all the elements in the first row of W 
arc T<nown. VJo can systematically make use of the relations 
(C.--.3), (C.4.5) - (6.4.9) and (6.2.20) to determine the 
oth.'r cloraonts of ¥. Let us illustrate this procedure by 
a simple oxamplo. 

rmo-^ric : For n=3, assuming that w^^, w^2 and w^^ 
arc ’;nov/n, BZpress other elements of W in terms of 
those throe. 

From (C.4,2) we can obtain 


w 

21 

= 1 

-f 

^11 

w 

= 1 



22 



12 

w 

= w 

13 


23 




(6.4.12) 
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» ^Zl * 

1 - w 

11 

^12 

35 - ^23 - 1 

'~^7' 


I 


w 


(6.4.13) 


iJsin:; (C.P.20 
b 


■Wp./w 
21^ 12 


'Z ~ -Wgg/Vj 


23 


(6.4.14) 


Combinific (6. ^.13) and (6.4.14' 




we can write that 


w, 


21 


w 


31 


( 1 -- w 


w ’ - "n ) 

12 12 


v; 


7 

3S 


^ * *82 - ’'si 



w 

w 

33 

= 

1 ■'^23 - 



^12 

^ • 4 * 

32) 

to compuue th< 


w 


13 


(6.4,15 


iCno'.vinr, thus the first two row elements of W, use 

-/ .4,15) tc compute the third row elements of W. Observe 

‘ r^.-^ (0.4.13), for b^ to be finite w =/ i, w ^0 

1 1 

ard •/ 0. ?or stable systems, such difficulties do 

not, arise. 
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6.5 ^CETHODS OF COWOTINO W FROM k 

The practical utility of the proposed canonical 
form depends upon the ease with which elements of W 
may "be determined from those of A.i At present the only 
use of this canonical form is in the construction of 
Liapunov functions, is the construction procedure 
through this canonical form any simpler than the usual 
method of solving the n(n+l)/S equations 

T 

A 9 ^ “ 2 ~ ~£ ? 

This is partly answered in this section "by presenting 
systematically a few methods of computing W. Also 
attempt has been made to find elements of W from the 
elements of the Jury table. This attempt is not 
completely successful. However, it may be possible in 
future to establish the link between the Jury table and 
the W matrix elements and to find some more uses for 
this canonical form as have been found for the Schwarz 
canonical form in the continuous case. 

Lot us consider first the straightforward method 
which is used while deriving W in section 6,2, First 
transform A to its companion form 0^. To obtain 

C from C , either we can apply scalar bilinear 

-d -a 

transformation, z = {s+l)/(s*l) to the characteristic 
polynomial F{z) of or we can apply the simple 

rules stated In [4?] which involve multiplication of 
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- jimod He trices. The latter approach is found 
no simple ond olefiant. However, to determine the 
ur u: .or 'tc.u ion matrix m [gee equation ( 6 . 2 . 12)7 , we 
ro n i. c to find r. Let the characteristic polynomial 


DO Liven cs: 


7.^ -f 


' ^^“1 + n „ n -2 

+ + . . . + a„ 1 z + a. 


(6.5,1) 


id ino 


1 + 

1 + a, + ar 

1 C. 


1+a, +ac + *»* +a^ 
12 n 


i} is biven ty |_45] : 


(6.5.E) 


P„ ( 5l “ ^ 

g n “1 


(6.5.3) 


\/hcro is a lower triangular matrix given by 


'ei’ii - V® 




■1 1 .1 


0 j > i 

P i > j 

n 


(6.5.4) 


ml D„ io given hy 


{P2)it 


Pn-j 


= 1,2, ... ,n) 


(6.5.5) 
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Having found p from (6.5.3), ^ can be obtained using 

the method discussed in fsoj , Note clearly that if the 
transformation S ( = P M N) is not required or if S 
is found by using the equation (6.3.22), then determination 
of D is not necessary, in such a case, directly 
can be found from explained earlier. Next form 

the Routh array of the characteristic polynomial of C^. , 
Use the relations between the bj^^s and the Routh array 
elements R..'s ['see chapter IlJ to compute the b.’s 
(i = 1,2,. .,,n). Apply the simple rules (discussed in 
section 6.3) to determine the elements of W knowing 
the b ^ ' s . 


Instead of using the rules given in section 6.3, 


recursive relations may also be used for computing Wj^^’s 
from These relations can be derived by simply 

equating the corresponding elements of (W B) and (W +'B h- I) 
see equation (6.4,1 )J , 

1 




w 


12 


(1 + Wjj) 


w 


22 


= 1 + w 


12 


(6.5.6) 


"^21 = -^^12 


To compute the elements on the first superdiagonal: 


w 


i(i+l) 




<'"ii " 


w 


i(i-l) 


), 


(1 = 2,3, ... ,n-l) 


(6.5.7) 
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I'D compute tJic elements on tJie second superdiagonal of W: 


w 


i(i+S) 


r — *1(1+1) - *ii] 

iti 


(6.5.8) 


• (i ** 1 } 2 j. , , , , n-2 ) 


To compute rest of the superaiagonal elements of W: 


w 


a 


t . 1 

1-1 


(i = 1,2, . ,n-3) 

(j = i+3, i+4,.,.,n) 


(6.5.9) 


To compute all tiie suMiagonal elements of W: 

(i-l) 

, j=l,2, , . . ,n-l, i=j+l, j+2 , , ,n 


w, . = w 
il li 



lc=j 


. (6.5.10) 


To compute all the diagonal elements of f; 
1 


w 


ii 


^i-1 • 


[*1(1-1) - '^1-1 - ’"l(l.s)],l=3,4,...,n (6.5.11) 


Let us also recall the common "ratio property j^For convenience, 
we call thus the property in (6.4.9) and (6.4.11) satisfied 
rospeotively by the superdiagonal elements and the 
suhdiagonal elements of discussed in section 6.4. 

It is true that the equations (6.5.6) - (6,5.10) 
do not explicitly toll how the computation of is ■ 

carried out from the knowledge of the hi’s. First let 
us assume that has been computed (we shall explain 
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how this is done a little later). We shall now careftilly 
explain how the other elements of W are computed. 

Find ‘^ 12 ,^ ^21 ^22 (6,5.6) and w^^ 

using (6.5,8) with i=l. Compute all the other elements 
in the first row using (6.5.9) i.e. 

1 


w 


l2 


— '" 1 ( 3 - 1 ) - '" 1 ( 3 - 2 ) ’ ^ 
3-1 


and then find all the elements in the first column using 
(6.5.10) i.e. 

^jl ^ ^li **^3-1^ » j=5,4,...,n (6,5.13) 

From now on there is a nice pattern of computation. Find 

using (6.5^7). As — ^ • = (common ratio 

23 Wj3 “in 

property), just knowing Wgg and the first row elements, 

determine all the elements to the right-of Wg^ in the 

second row. Use (6.5.10) to find Wgg. i.e, W32 “ *'^2 ‘'''^23* 


w 


32 ^42 


As 


^!ii2 


w 


31 


w 


41 


w 


nl 


second column helow w^ 


, find all: the elements in the 
Now find w „ from (6, 5 *11). 

32 00 

Continue this process of computation exactly in a similar 
pattern imtil the element w^ is computed. 

The above described computational scheme is 
recursive in nature and this is well-suited for use on 
computer. Let us next see how we can compute Lrom 


W = I + 2(B-i)‘ 



]58 


we obtain 

w « 1 + 
11 


2 [^c of 86 tor of 11 element of (B - I)] 


L Determinant of (B - iQ 


It is not very difficult to see (Refer Appendix I) that 


« 1 + 2B(n-l)/B(n) (6.5.14) 

Now B(n) and B(n-l) can be easily obtained from 
the recursive relations (Appendix I) 


B(o) = 1 


B(l) 

B(k) 


-(1 + 

-B(lc-l) + B(k-8) 

k = 2 ,3, . . . ,n 



(6.5.15) 


So we infer that the rules for determining from 

bjj^’s discussed in section 6.3 are useful for hand 
calculation whereas the above described recursive 
computational scheme is well-suited for use on computer. 

In the rest of this section, attempt is made to 
relate the Schwarz elements b^^’s with the elements of 
the matrix A which is assumed to be in the phase- variable 
form and to relate the elements of W , , with those 

of i.m It is also speculated that such relations may bo 
expressed through the elements of the Jury table of the 
characteristic polynomial of A. Even thoiagh complete 
success is not obtained in this investigation, it is hoped 
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that the results to be presented ■will be of considerable 
use for further investigation. Before proceeding further 
let us recapitulate that in the continuous case, the 
.elements of the Schwarz matrix are related to the elements 
of the characteristic polynomial (or the elements of the 
corresponding -phase-variable form) through the Hurwitz 
determinants or through the first column elejpents of the 
Routh array. Thus intuitively it is felt that the 
elements of W matrix may be related to those of A 
matrix through Schur-Cohn determinants or Jury table 
elements. 

Let the characteristic equation of A be given 
by (6.5.1). Let J. . denote the ij-th element of the 

X J 

Jury table of (6.5.1). The following results are 



(6.5.16) 


(6.5.17) 
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n 4 

(4 4' 2a 2 - Bag - 434 . ) 

*1^ sg n il, Iiiiii mi ii,ii i. ' 

4 . U + + ag 4 . + a^) 

(20 - 20 a| 4 lOa^ - 18ag - ISa^ - 2 a^ + Bagag 
- lOaga^ + lEa^a^ + ISa^a^ + 2a^ - 2a^a^) 



(1 4- ag + ^g + ^4^4 + 2a^ - Za^ “ ^®4^ 

(64 “ 64a, — 64a_ *■ 64a^ + 128a_a + 64a a_ 
4 ^ o 4: 1 o 

+ 64a,, a„a. - 64afa. - 64a. + 64a^ - 64a_a^) 
1 3 4 1 4 4 4 2 4 




(20 - 

20 a| + 

lOa^ 

- 18 ag 

*r 

-12 ag 

4 + 2a^ - 

sag - 4 a^) 

- 2 a§ 

+ 23383 

- 103334 + 

128284 



^ *4* 2B^ - 

Bajag) 

- 

(1 - a^ + 

ag - ^3 + ^ 

4 ) 

** 

Bag 

- 


(20 - 20a| + lOa^ - 

ISa^- 

IBag - 

0 2 



- 10 V 4 

+ 12a 2^4 + 

ISaTa 

1 

4 + 

- 2a 

1®2* , 







• • • (6 

.5.18) 


For liiglier order systejns the corresponding expressions become 
more and more oomplioated. In the above, are related 

to 84 ’s directly. We shall now express these relations tbrotgh 
the Jm-y table elements. We shall denote the 
characteristic polynomial of a by F(z). ' 



n = 2 


161 


2 (J-, 


21 “ 


b, = / T \2 



21 “ ^1 J + (^ 



" - ^3i)J 


F(l) [ 3 (J, 


31 ■ -^ll) ■* (Jss - JjgO 






(6.5,19) 


(6.5.20) 


bj - (-1)= SldJ_ _1 

F(l) 

-'■or i:.4, n=5 eto. , tHe expressions for b ana b 
Hhow so» pattern wbereas for otbers. it is not possLle 
present to laentlfy any suob pattern. Expressions for 
*’n bj for n=4 ana n=5 are giren below; 

n = 4 


^‘^21 ■ ^11> ■" 2(^28 - ^ig) 


(- 1 ) 


r(i) 

4 linii _1_ 


> 

»«/ " 

1 

n = 5 


I’d) 


(6.5.21) 




5(J 


31 ■ * (J'p.s - ^,,1 


23 "13' 


F(l) 


:i = (-1)5 gJ,~l) __L- 

^ F(l) bgbg 


(6.5.22) 





Let us now 
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i i'll simple cases. 


n = 3 


■ h * 


1 ~ l)/3 

-(1 + a 

+ ag)/8 

(1 - a 


1 * ®2 )/2 
1 “ ^2 ^ 


(6.5.23' 


n = 3 


'3 - 2(1 - ag - a2 + a^a^) 


(3 + a 


2 - Sg - 5a ) 


- a .)/2 -.1 


“(1 + a + a + a_)/4 


1 2 °3 
(®3 “ ®2^ ^ 


(1 *. Bo - + a ) 


2 (3 + - ag - 3a ) 


^^3 “ 


- ( 1 +• a^ + Sg + a2)/4 




(Eg - a ^) 


2(1 - a - a^ + a.a_) 


[6.5.24) 


h ^ ^ ' 2 °5 " ^1°3^ ^ 

2 (3 + aj - ag - Sa^) J 

( - ^2 - ^3 ^ ^^3^ 

l.(-l + (a^ - aj)/2) (3 + - Sg - Sa^) 

2(1 - ag - + a ^ ag ) 

(3 4 8 | - Sg - 2 ^ 3 ) 

3(1 + - a^ - a^) 

3 * — 

(3 + ai - ag - 3ag) 



163 


may be recalled (section ^ t 
of w that only n elemAnt 

^ are independent q, elements 

to be tbo ■ ‘ “ 

° De those in the f4-n«4- 

i^inst row of w 'nhan tu 

iollowin^ reesnu- „ “* fhon the 

^ins results are obtained 

a = S - 

r(i) 


W 


18 


w 


n = g 
w 

13 


w 


18 


’ll 


- ill) 


^11 “ ^ 

, 88 


^il - ^83 


( 6 , 5.85 


w 


1-1 




" ( 6 . 5.86 


Even though such results are not 

Mgber oraer svet„ P^eeeatly amiable for 

y ms, we obseyTe that even In hlahev 

°raer systems 


w 


In “ - 



( 6 . 5 . 87 ) 


In what follows, an attempt Is maae to o 
Seaeral the first r , 

elem t tlie 

slements a* tg . 

, -■ ®uoh an erpresslon 

involves some of the Bonth • 

uth array elements of the bilinear 
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ir -'nr.iurMed nystepi. However, tlie expression obtained 
is J’lioc iiiueresGing and may be of use for further 
i nvcstication. 

The mobrix a is assumed to be in pha-ge-variable 
''o 'm. Tlio expression to be given is obtained after 
'^ener'li v ng the results from systems of order upto 
sc“'on. "he method given in the beginning of this section 
nsed Pox obtaining the results for low order systems. 
Bef'oro stiting the results, let us define the follovsring 
ou TtJ tins . 

r 

f ^ ^ (6.5.28) 


Wlf'' 

Ij 


ij-j nxn 

,j , j = ],2,..,,n 

^ , i > 0 

+ 2F, 


0 


0 


iJ 

^Tj - '■'kJ-D ®hi-l)(3-l)’ 

liitv th'<t F is an upper triangular matrix. 


0 0 ... 0 Pj^/2 


( n+1 ) xn 


,'!if ro 


P, 


1 + + ag + a^ + • ’ *+ a 


n 


(6.5.29) 


(6.5.30) 
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•i+1) 


' > - = 1,2, ...,n 

> J -^ ( n-i+i ) , i > 1 


(6. 5. SI) 


i( j+1) 


+ S 


(i-l)(3+l) 






(i .0 
v.'lth 


tio elements of s in a tow alternate starting 
the croos-diagonol entries and moving leftwards in 
row. rote that s is a cross-diagonal matrix 
. rlo-r.ents lying below the cross-diagonal are zeros) 
-tl the cross-diagonal elements equal to -2. 


(6.5.32) 


0 J nx(n+l) 


1 T '1 

(n-2k)x(n-2k) 


n - 2k 


I (w-i' 


1 , ' = 1 ,2, ... ,m 

C , j > n-i+1, i >1 


■1(J^1 


(i-l)(j^-l 


+ m-j+1), i= ift-l, m-2,...,l 


i + T 




, J <m-i+l, i < (m-1) I 

y 

... (6.5. 33 ) 



Oi 
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n - the elements of in a row alternate starting 
i rvT. I he cross-diagonal elements and moving leftwards in 
ti. .t rev,. 'lote that also is a cross-diagonal matrix 

• M,'t thi, cross-diagonal elements equal to 1. 


n-2k 


-2k 




?k 


0 


8k+l 

0 


0 


nx (n+1) 


(6.5.34) 


fi 


‘n-1 


’n-2 


1 , 1 ) 



n-2k 


n-2k 


2k 

0 


2k 



0 


nxn 


(6.5.35) 


(6,5.36) 


v/hrrn p io an (n-2k) x (n-2k) diagonal matrix, diagonal 

-K 

onlrior lelng determined by the Routh array elements of 

^ ( one ocction 6 .2 ) . 

’Vith the above nuantities, we are now ready to 

prc.k.ut tho relationship between iw and a. 


w 


in-1 ) 


+ 2 



F a 


(6.5.37) 


•whui'i- = 2 

-0 

centuinod in k . 


and 


[^k"] means the largest integral 


The equation (6.6.37) may be viewed as 



iCii 
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the transformation between w and a. 

1, viewpoint becomes more apparent if we rewrite (6.5.37) 


3 



't' .ref’Oly that the transformation matrix K involves 
/ 1 ^ 


(l - 0,1,..., 


ij ) and I , that these matrices 

2 

.j'fj r*, rmva <r?itc simply and that they contain Just numbers 
i.f . ihoy ao not depend upon the system. The dependence 
n tic- system comes through in (6,5.38). 

Let us illustrate these ideas by a simple 


, . pqp n = 4, determine K. 

1 A A A ^ 

y = i (S + 2R, T ) F 
- 0 -1 — 

Ihf relation w = K a. is illustrated in the 
ncjxl i -00. 

Blurents of in (6.5.36) may also be expressed 

**jcC 


In -onorvl. 
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6.6 HiOBI^.rOS fOTTRS IMMSTIOATIOII 

*^®“^^^,f,#Sousaiiig ia tJie:.ppeYlons sections 
several woWe,M for farther lavestigstion asVe been posed' 

la tais^seOtlon. ,^^ltlon to presenting' th^^p^oble^ ‘ 

'/ V/, 4;-' " '''* t of estabiigJbiijQg a iiak *' 

gtabllltv oflte-rta #aa /tnl a m . 

c-rxTjeria X48 , 49j ©nd Liapunov’s 

*’® °oasl6ered. 

0 “"y Be given for ensuring 

r. that the.negatlve saMdeflmteness of Q ip implies 

the asymptotic stability of the system In (s.sip). Second. 

the relations between the elements Of the W matriz and 

those Of the a matrix may be established In general 
either through the S=hur-Oohn determinants or through the 
elements of the Jury table. Third, some more uses for 
this oanonioal matrix w may be found. For example, 
analogous to the Schwarz matrix being used for performance 
measure evaluation In the continuous case, w also may 
be put to such t use. 

Finally, let us consider the problem of establishing 
a link between the Jury’s stability constraints and 
Liapunov’s second method, m section, 6.3, it has been 
found [see equations (6.3.25) - (6.3.27)] that the^ -i. 
application of the second method of Liapunov requires,"'^* 
for the asymptotic stability of the system, the following, 
conditions to be satisfied by the elements of W: 
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0 » k = 


( 6 . 6 . 1 ) 


1 w + w 

^(n-1) nn 

1 - w 


(6.6.S) 


jjGt 'ir. oLo how these constraints can be shown to be 
n".i.i\f?jlsnt to the Jiiry’s stability constraints 1^4^ , 
i't-r tho cases v;hen n=2 and. n=3. 
n ■■ 

/rom (6,5,83) it is not very difficult to see that 


end a. 


-(*11 -^* 28 ' I 
<*11 *88 - *18 *81 > 


(6.6.3) 


aIsc Irom 


B + V/ + I 


VC rr-*i (coe section 6.4) 


w^i = 1 f w^^ 

V/gg = 1 + w^g 


(6.6.4) 


’ury’s stability constraints for this problem are: 

.’{!) ) 0 , F(-l) )■ 0 and 6 0 (6.6,5) 


?(1) >0 -^^l + a^4'ag> 0 


(6.6.6) 


(C.6.3), (6 


.6.4) ^nd (6.6.6), after some manipulations 






equiyalent 


, - , J'*' ^ ■*■ Til 

la view Of (6.6.4) this aeata 


-"' i 2/*21 > 0 


Owing to (6.6. 5 J an! (5.6,4) 


Combining ( 6 . 6 . 9 ) and (6.6!u) 


this, in view of (6.6.10), implies 


lubstmting (6.6.2) and (e.a^i) in (e.eilE) 
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Obtain 


,(6^6.14) are_ e;sactl 
jj putting n*2 in (6.6.1) 


ll*- e^aiiDXe,, we started with 
constraints and then after many maj 
tiiat tbese are equivalent to tbe sl 
s obtained by Liapunov’s second met 
ceed in tbe reverse way, por n=3 


Ifelng the results 
(6.6.15) 


can be. shown that 


(6,6.17) Imply 
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■ ‘ ^ ' fi:0l 

J’(l) > 0 


?(- l ) < 0 


' -J, ■■■.■., . 

1 - ag ^ 0 


( 6.6 


ilury’s constraints [487 are exaotUr ^ -i ^ 

L oj are exactly similar to those 

i-H (6.6.37) except for the far»+ ■(-hQ+. u 

the fact that he gets in addltioi 

to these ' 

^ ® 3 0 

(Actually he gives a condition )a^| ■< i., ^ 1 ^ 
would Imply 1 . a,j > 0 , 1 . > 0 ) . Dees this mean 

the condition 1 + a^ > q is not required ? Of 
course, this needs further Investigation. We end this 
section hy saying that it ®y be possible in general to 
establish a link between Jury's stability .constraints 
and the second method of Liapnnov. ’ V" ' 


27) 



Ca^rTER - VII 

COVCLESIOIE 
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The canonical form introduced by Schwarz serves 

t 

Os a useful way of characterizing the performance of 
linear tlme-invarianb continuous-time systems. It has 
several useful properties which have been listed in 
this report. It is also shown to be useful in 
simplifying large dynamic systems and in deriving a 
canonical form for linear time-invariant discrete-time 
systems. The problem of transforming any given system 
to this canonical form has been dealt with very 
systematically, a Vandermonde-like matrix which 
transforms the Schwarz form to the Jordan form has also 
been derived. The Schwarz canonical form has been 
extended to multivariable systems also. A simple network 
interpretation has been given to the Schw'arz canonical 
form. 

The method for simplifying large dynamic systems 
investigated in this thesis is a matrix generalization 
of the associated transfer function introduced by 
Gustafson . However, the approximation made is 

very much influenced by the type of inputs applied to 
the system. Tt may be possible to devise a method for 
simplifying large systems, which is optimal with respect 
to some particular input. This may call for the application 



this direction has been recently made by Meier and 
Luenberser I_38j . 

A new canonical form for linear time-invariant 
discrete- time systems introduced in this investigation 
has many useful properties similar to the Schwarz form 
in the continuous case. Hovrever, the investigation 
remains somewhat inconclusive s^nce the relations between 
the canonical matrix elements and the Jury table elements 
are not generalized. If this is done, the Jury table test 
for stability may be proved through the second method of 
Liapunov. 

Throughout the presentation of the results in 
this thesis, wherever necessary, problems which require 
further investigation have been posed. A brief summary 
of these problems will be given here, 

1. A technique of modeling the physical system may be 
devised to yield the system description directly in 
the Schwarz form. 

2. The Schwarz canonical form for multivariable systems 
may be investigated further and its practical uses 
may be determined. 

3. Tor the simplification procedure introduced in this 
thesis, estimates of error between the actual and the 
approximate S 3 ^stcm responses may be determined. Some 
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difficulties associated with, systems with complex 
eigenvalues maj be resolved. It may be possible 
to find out a performance measure which becomes 
optimum for this simplification procedure. The 
problem of extending this method to multivariable 
syst'-ims also needs further investigation. 

4.. In connection with the proposed canonical form 
for linear time -invariant discrete-time systems, 
the problem of relating its elements to Jury 
table needs further investigation. 

5. In the case of interconnected systems, the 

investigation of the stability of the entire system 
through that of subsystems has been carried out 
recently [50j , It is interesting to ask whether 
there? will be any simplification if the individual 
subsystems are in Schwarz form. It may be possible 
to construct a Liapunov function for the entire 
system through the Liapunov functions of individual 
subsystems at least for simple modes of 
interconnection pattern. 



^PSMDIX - I 

INVERSE OF TRIDlAGOmi MATRICES 
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Let A be a tridiagonal raatrix. That is, if 


then 



A 



ran 


- 0 when D 7 ^ i or 2 ^ i-1 or j i + 1 


Here, it is aimed to state the general expressions for 
the elements of A-1 in terms of explicitly. 

These expressions are obtained by generalizing the results 
got from low order examples. Let 


L xjj jixn 

The resirlts will be stated now fop three cases. 


1. 

When 

A 

is 

any tridiagonal matrix 

2. 

When 

A 

= B 

where B is the Schwarz matrix 

3. 

When 

A 

= (B 

- I) 


Case 1; The elements 0..’s are given by 

^ t) 


li B{n) 


F(i-l) B(n-i) , i = 1,2, ...,n 


i-1 


^ F(j-l) B(n-i) (-1)^-*^ > 3 


C. . 

il 


B(n] 


F(i-l) B(n-j) (-l)^*^^' 


Lili 

r\.®k(k+l) 

k=i 




» j > i 



where 


B{0) 

E(l) 

B(k) 


= 1 


= a 


and 

F(0) 

F(l) 

F(k) 


nn 


’(n-k+l)(n-k-!-l) ®(k-l) “ ^n-k-^S) (n-k+1)^ 

® ( ii-k-*-2 ) ( n-k+3 ) ® ^ ^ 

k = 2,3j,,,jn 


11 

\k ^ 


k(k-i) ^k-i)k 


k = 2,3, , ,n 


.. 1 . 1 ) 


Note that B(k) is the k-th ascending principal 
minor and F(k) is the k-th descending principal 
minor. Clearly 

F(n) = B(n) 


Case 2 ; 

C 


When A = B 


ii 


f 0 i even 

1 — 1 — p(i_i) B(n-i), i odd 


LB(n) 


1 — 
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0 


J even 
(i-1) 


1 j 

i-/' 


J odd 


0 


even 

bT^ r(l-l) B(n-J) (-i)iV ) 

i Odd 

• •• jH 


Oi 


'viiere 


Hi) = F(k-E), ic = 2,3,...,n_i 

F(0) = 1 


F(l) = 0 
and B(k) = b 

B(0) = 


n-k+1 ®(k-2), k - 2 , 3 j,..j 


n 


B(l) = -b 


n 


Note iiere also that F(n) = B(n) 


• { » 1 • 2 ) 


Case 5 : 

= (B - I) 


^ii= 


Bin) 


F(i-l) E(n-i) , i = 1 , 2 ,..., 


n 


B(n) 


C. = 


Bin) 


^ B(i-i) B(n-i) , i > 

tk=j J 


F(i-l) B(n-j) j 


i 1, 1,J = 1,2,..., 


n 


vjhore 
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F{0) = 1 

F(l) = -1 

I’(k) = -F(k-l) + F(k-S) 

k = 2,3,...,(n_i)jii 

and 

B(o) = 1 

B(l) = -(1 + b^) 

B(k) = -B(k-l) + B(k-2) 

k = 2,3,...,(n-l),n (ji.1.3) 

Note that B(n) = F(n) , 


APPllIDIX - II 
MDRM^LIZATION PROCEDURES 

Two normalization procediires are discussed here, 

I. Consider the transfer function 


G'(s) = 


z:(s) 

v(s) 


(A.2.1) 


Cos“ * 


+ • 


+ C^_^s + 


Divide the denominator and numerator of the expression 
in the right hand ©ide of (a. 2,1) by Cq. Define 


&/C, 


(A. 2,2) 


i — lj2j,,,jn 


0 



We get 
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G^(s) = 2l s , ). 

v(s) 


+ a 3 ^“^+...+ a„ ,s + a 
1 n-l’^ n 


(A. 2. 3) 


Define 


Then 


bv( s) 


u(s) 


(A. 2 * 4 ) 


- sill 


^ *...* 


(A.2.5) 


G(s) is tile normalized transfer function used in section 

5.2. 


II. Consider the differential equation 


+ ^ S 4'...+ C , ^ + O x = g u(t} (A. 2. 6) 


Substituting X = CO t where do = n J ~ iu (a, 2, 6) 

and dividing the resulting equation both sides by 

C • 

and replacing ^ CO^~^ by (i = l,2,...,n-l) and 

-[jy y{t) yield 
n 


SIe + a, Alii- * . . . + ^ = v(t) 

d^n 


(A.2*7) 


Taking Laplace transforms on both sides of (A. 2. 7) 


3c(s) 

^(s) 


(4.2,8) 


s^ + a s^-^ +••*+ s 1 



'rhis is tiae normalized transfer ftmction corresponding 

(.^.a.6). 
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APPENDIi: - III 
NIBffiRATOR DYNA^^IC SYSTM^ 

It is desired to extend the method of simplifying 
large systems, proposed in chapter Y, to systems with 
numerator dynamics. Consider the system given by 

x^^^ + +•••+ x^^^ + a^^x 

= d -fdn d u (A#3,l) 

0 

where m can be at most eoual to -1 , where ^ is the 
size of the simplified system. (If we were to work with 
the system (a. 3.1) Itself without worrying about any 
simplification, then m can be less than or equal to n) . 

Two methods will be given to reduce (a. 3,1) to the 
required form (5,2.3). These methods are discussed in 
detail in q43 [ , 

Method 1 ; 

The system given in (A. 3.1) can be equivalently 
represented in a matrix form as under: 
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y- 


^3 


0 

0 

0 


1 0 
0 1 
0 0 


[ynl 


-a 


where 


n 


-a 


n-1 “2 


n-2 


X - Oo u 


... 0 

0 


P^' 


'h' 

• • • 0 

0 


^2 



• 

• 

♦ 

o 

0 


^3 



• « • • 

« 

j 

« 

+ 

• 

* • # * 

9 


• 


• 

• • • • 

• 


• 


• 

. .. -ag 

-uj 

L 



-^n,. 


u 


(A. 3. 5) 


am if a = n , then 


P° 


1 

0 

0 

0 

• ♦ • 0 

c 



®1 

1 

0 

0 

• . . 0 

0 



^2 

®1 

> 1 

0 

• • • 0 

1 

Co 



®3 

^2 

"l 

1 

... 0 

^3 

• 


• 

• 

• 

• 

• ♦ • • 

• 

«t 


• 

• 

• 

• 

• • • • 

• 

■ { 


a.-, 

n 

^n-1 

n~2 

S-3 

* • • X 

^n 






J J 

- J 


If m = -ll 4 n 


(A,3.4) 


then 


= di = dp = d, 

mu 0 12 

Thus knowing dj,'s (i = 0,1,..., n) 


= 0 . 


and 


ai’s (i = 1,2, ...,n) 


can find c^’s by inverting the (n+l)z(n+l) matrix 
in (A.3.4) and premultiplying the vector in the left hand 
of (A. 3, 4). Again another coordinate transformation 
can be applied to take the system description from the 
form in (a. 3.2) to that in (5.2.3). Note also that the 
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initial conditions for the vector i = (y 3 _ ,yg ,...,y^)‘^ 
in (A. 3, 2) are given by: 


o 


'x(o) 


^0 

0 0 

.. . 0 


u{o) 

73(0) 


x(o) 



0 

0 

0 

.. . 0 


ft(o) 

73(0) 


x(o) 



0 

M 

0 

0 

... 0 


u(o) 

• 


• 


m 

« • 

* • • • 


« 



1 * 


• 

• • 

« • • • 


• 

« 


• 


* 

• » 

4 • • * 


• 

• 


• 


1 • 

* • 

• « • » 



0 

1 



i 

Vi 

‘^n-S °n -3 

... Cq 


„(n-ir 

" (o)_ 


Method II: 


By this method, the system in (A, 3,1} can be 

equivalently represented in the form given in (5,3,3) 

where now the relation between x and x, , Xo , ... ,x 

1 ^ n 

is given by: 

X = b^u (dj^ - ajoid^) x^ + (*^11-1 “®n-l‘^o^ ^2 '*'*** 

.... (dg -agd^) . (d^ (-1.3.6)* 

To determine the initial conditions for (5.2.3), 
differentiate (a. 3. 6) success^ '"^ly (n-1) times each time 
substituting for from (5.3.3), At t=tQ 


'•'^Note that in writing down (^.3.6) » 
that ’m’ aopearing in (ji,.3,l) equals n. 


then d. 


= do = 


= d. 


n- 


= 0 . 


it is assumed 
If m = n, 


• w » 
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(initial time), (A. 3, 6) and the (n-1) differentiated 
equations v;ill give the necessary n equations to 
evaluate . ,zjt J . 


n' 


As far as the simplification taethod suggested 
in this work is concerned, method li seems to be easier 
than method I, in section 5.5 a simple example system with 
numerator dynamics with n = 1 is reduced to a system 
of order two. 


APPENDIX - IV 

SIICPLIEICATION OE A SECOND ORDER SYSTET!d 

In section 2 of chapter V, it is stated that 

unless / bj^, the comparison process must be 

performed till b is compared with b . Suppose the 

^ 1 

comparison process has ended with bg being much larger 
than b^. We infer that the given system may be reduced 
to a second order system. Can we set up a criterion to 
check whether this reduced second order sysuem can be 
simplified any further (i.e. can it be reduced to a 
first order system) ’ Obviously no one likes to reduce 
a very high order system to a first order systemi Yet 
for completeness sake we v;ill answer the above question. 

It is very obvious that the second order system to 
be reduced to a first order system cannot have repeated cr 
complex eigenvalues. Let us now state the required 
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criterion. Whenever "*^2 

r 

b^ £. bg but each of them being much larger than 1, 
infer that the given seconi order system can be reduced 
to a first order system. 

it 

Sr ample s . 

1) Given a system as in (5.2.3) with a^^ = 101 and 
Hg = 100, Therefore, b^ = 101 and ~ 1^0. 

Here, = "bg "^2 this 

system can be reduced to a first order system. In 

% 

fact the eigenvalues of the given system are -1 
and -100. 

2) ^2 = 3^ = 20.1; b^ = ag = 2. Here bg ;^> b^, 

In fact the eigenvalues of the system are -0,1 
and -20, 

3) bg = ai = 1010; b^ = 10,000. Here bj_ » bg. 

The eigenvalues of the system are -10 and -1000. 

4) b = a, = 3; b^ = a„ = 2, Here b, = bg and 

g i 1 li 

b^ and bg are close to 1. Therefore no 
approximation can be made, (The eigenvalues of the 
original system are -1 and -2), 

5) If the system has repeated or complex eigenvalues, 

this criterion fails, 

(i) tg = = 200; ^^=83= 10,000. Here 

b >> bo. But yet eigenvalues of the system 
1 2 

dxo -100 8iid -100# 
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(ii) bg = = 1; ^1 = ^2 = 25.25. Here "^S* 

However, the eigenvalues of tij-is system are 

-0.5 + j 5. 


So in conclusion it may he stated that whenever 

we have hg = h^ and both bg and b^ are much 

larger than one, reduction can be performed. The 

reduced first order system has b^ = — . in other 

1 b2 

cases, it is nece ssary to find the eigenvalues 

bo + x/bo -4b 1 n 

-H_= ^ i. I and if they are sufficiently 

2 

apart, oerform the reduction. 




APPEl^IX - V 

ELIMimilOH OF REDIM)ANT STATE TARIaBLES THROIXIH THE USE 
OF SCHMZ CaNONIC.IL FORM 

When several subsystems are assembled to form 
the main system, it is quite possible for such an 
assembled system to have somo redundant state variables 
[36, 44l . Presence of such redundant state variables 

can bo easily detected when the main system equations are 
transformed to the Schwarz form. If the system is stable, 
then zero values of b. ’s would indicate the presence of 
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T edundant state variables, if there is one redundancy 
then b^ will be zero, if tvjo then both b-, and b„ 
will be zero, and so on. This state'Hent can be very 
easily proved. Redundant state variables introduce zero 
eigenvalues. Presence of zero eigenvalues can be detected 
fron zero values of bottom elements of the Routh array. 
That is, presence of one zero eigenvalue, for instance, 
results in =0, Expressing b^^s in terms 

of the Routh array elements 

h = 

^ ^(n-l)l 


^(n-2)l 


Hence the proof. 

The redundant state variables can be easily 
eliminated by discarding the first ro^<s and columns of 
the Schwarz matrix and also of the transformation matrix 
P. 7fe shall illustrate this by means of a simple example. 
Consider a system 

^1 
Xo 


X. 


T 

“41 


\ 

o 

o 

o 

i 


"^l' 


'o ' 

0 0 10 


Xo 


0 




4 


0 0 0 1 


^3 


0 

-®4 -Qg -^2 


^4 


1 


a - 




u 


i.C, X = n X + 


u 



1 


0 


0 


p 



1 

0 


bi +b 
1 2 


0 

1 

0 


0 

0 

0 


Let the transformea 


system U given by 



Suppose = 0. 
^‘irst row of b 
■ to be eliminated 



Then aieoard tfe eolumn ana 

ana the oorrespoMi^ Tariable z 
• This Is lllustKtea below. ^ 


iias 


-Zj- 


1 Q q’ 


^2 


o 

r-t 

0 

1 


^3 


j> ~bg 0 1 




-P “ -^3 -L, 



The transformation matrix p 
its first row and first 

— i 

-4 Q 0 Q 


.Bf 


p. ^ 

-0- 

^8 


' 0 

^3 


0 

^4_ 


1 





is also to be changed 
oolunmhave to be di 


similarly 

scarded. 


P = 



6 


1 

0 


b-, +b,. 

, 1 2 

^ =bp 


0 

1 

0 


0 

0 

1 



191 


LIST 01' REPEInCES 


1. H.R. Schwarz, ”A method lop heterminirg stability 
of matrix differential equations", Z. Angrew. Math. 
Phys. (in German), Vol.7,jp ^3-500, 1956, 

2. R.E. Kalman and J.E. Bertijn^ "Control system 
analysis and design via the second method of 
Liapunov,” Trans, ASME, Ser. d, pp 371-393, June 1960. 

3. p.C. parks, "A new proof ol the Routh-Hurwitz 
stability criterion using the second method of 
Liapunov”, Proc. Cambridge ?hii. ‘Soc. , Yol,58, 
part 4, pp 694-702, October 1962. 

4. P.C. parks, ”A new'proof of the Hurwitz stability 
criterion by the second method of Liapunov with 
applications to "optimum” transfer functions”, 

Joint Automatic Control Conf,, pp 471-476, 1963, 

5. J.l, Diamessis, new method for calculating system 
performance measures”, TEES Proc. (Correspondence), 
Yol.52, p 1240, October 1964, 

6. C.E. Chen and J.E. Diamessis, "Obtaining performance 
measures from Routh’s algorithm", IEEE Proc. 
(Correspondence), Vol, 53, pp 659-660, June 1965. 

7. S. Barnett and C. Storey, "Solution of the Liapunov 
matrix equation”, Electronics Letters, Yol.2, pp 466- 
467, December 1966. 

8. H.M. Power, "Solution of Liapunov matrix equation 
for continuous systems via Schwarz and Routh 
Canonical forms”. Electronics Letters, Yol,3, 

pp 81-82, February 1967. 

9. H.M. Power, "Further comments on the Liapunov 
matrix equations”, Electronics Letter , Yol*3, 
pp 153-154, April 1967 



^ 'I Liapunov matrix 

a^diagon^ input matrix, obtained 
^«trix inversioiH Electronics Letters, 

Vol.3, pp 325-326, Julyi967. 

Stores »Tiie Liapunov matrix 
equation and Soiiwarz's IEEE Trans, on 

Automatic Control (Correspondence), Vol. AC-IS, 

PP 117-118, February 19< 

1C. 3. Barnett and C. Storey, ?«sonie applications of 

the Liapunov matrix equation", J. Inst. Maths, 
applies, Yol.4, pp 33-42,1968. 

13. G.F. Chen and h, Chu, 'Li. matrix for evaluating 
Schwarz's form”, IEEE Traas. on Automatic Control 
(Short Papers), Yol ..iG-ll, pp 303-305, April 1966. 

14. C.F.^Ohen and H. Chu, 'Li. general proof of the 

matrix and its inverse", ibee Trans, on ^ 

.automatic Control (Correspondence), Yol.i-iC-lS, 
pp 458-461, august 1967. 

15. D.G.^Luenherger, "Canonical forms for linear 
multivariable systems”, IEEE Trans, on automatic 
Control (Short Papers), Yol. aC-lB, pp 290-S93, 

June 1967. 

16. H.S. Wall, Analytical theory of continued fractions 
(Book),D. Yon Nostrand, New- York, N.Y. , 1948, pages 65, 

173-176 and 182. 

17. J. Maofarlane, "The calculation of functionals 
of the time and frequency response of a linear 
constant coefficient dynamical system'’ , Qu^ft. J. 

Mech. Appl. Math., Yol,14, pp S59-271, 1963, 

18. P.C. Parks, "Stability analysis for linear and 
nonlinear systems, using Liapunov’s second method", 
Progress in Control Engineering (Book), Yol,2,* 
Heywood & Company Ltd., London, pp 31-64, 1964, 

19. IEEE Trans, on ..o-utomatic Control (Correspondence), 

Yol. aC- 11, pp 607-610, July 1966. 

20. H.Ti:. Power, "The companion matriz and Liapunov 
functions for linear multivariable time-invariant 
systems", J. Franklin institute, Yol, 283, pp 214- 
234, March 1967, 



21 . 


193 


22 , 

83. 


24. 

25. 

26 . 


27. 

28 . 

29. 

30. 

31. 


32. 

33. 


teaJsfSmat-l n; ; Tlswanathan, "On the 

DubllehJf fi fa °x®°***'* canonical form", to be 
SutniSf? tile June 1 \q issue of ISEE Trans, on 
automatic Control (Corj>gpondence) 

varlah?P^ tansformation to (phase- 

OnntSnT' r ^ Trans. on Automatic 

ontrol ( Correspondence ),7ol .AC-11, p 607, July 1966. 

Tiswanatiian, "a note 
fi-oi Schwarz to Jordan 
» ‘fco P^tlished in the June 1968 
issue or IEEE Trans, on ^^tomatic Control (Correspondence) 


’’Determination of%e inverse Vandermonde 
matrix”, IEEE Trans, on «.ntomatic Control 
(Correspondence), Vol.i^C-9,p 314 , July 1964. 


Custafson, ’hi paper anci pencil control system 
design , Trans. .iSME, Ser.D.^ pp 329-336, June 1966. 


M.E. van Valkenburg, Introduction to modern network 
synthesis (Book), John Vfiley & Sons Inc., New York, 
pp 397-400, 1960. 


T.G, Jdarshall, ^Ladder networis synthesis by 
similarity transformation" , preprint , Proo. w-llerton 
Conference on Circuit and System Theory, October 1965. 

C,H. Wilts, Principles of feedback control (Book), 
Addison Wesley, Reading,' Hass., p 66, 1960. 

Israel Navot, ’’The synthesis of certain subclasses 
of tridiagonal matrices with prescribed eigenvalues”, 
SIaM J. .ippl. Hath., Vol.15, pp 241-251, March 1967, 

T.G. Marshall, "Primitive matrices for doubly 
terminated ladder networks”, Preprint, Proo, nllerton 
Conference on Circuit and System Theory , October 1966. 

T.G. ^Marshall, "Tridiagonalization and related 
matrix and network algorithms for a polynomial space”. 
Preprint, Proc, itllerton Conference on Circuit and 
System Theory, October 1967. 

J.G. Truxal, Control Engineer's Handbook (Book), 
McGraw-Hill Book Co,,lnc., p E-14, 1958. 

E.J. Davison, ”4». method for slmp^fying linear 
dynamic systsms”, IEEE Trans, on Automatic Control, 
Vol. ^iC-11, pp 93-101, January 1966, 



194 


34. 


35, 


36. 


37. 


38. 


3 '. 


40. 


41. 

42. 

43. 

44. 

45. 


46 . 


I’.R. Ohidambara and E.G. Davison, "On a meV 
far simplifying linear dynamic systems”, IlElv 
T'ans,' on Automatic Control (Correspondence), 
V(1 .a(3-12, pp ii9-l£l, February 1967. \ 


M.l, Chidambara and E.J, Davison, ‘'Further remarks 

on simplifying linear dynamic systems”, IEEE Drans. 
on 'Automatic Control (Gorresnondence) , Tol.AC-lS, 
pp U3-S14, April 1967. 


R.D. Gustafson, 'Simplification of state-variable 
systms for rapid digital simulation", submitted to 
the 968 JACC. 


M, Acci, "Control of large dynamic systems by 
aggregations", to appear in IEEE Trans, on Automatic 
fontri, June 1968. 

Lewis \|e5er III and D.G. Luenberger, "Approximation 
of linar constant systems", IEEE Trans. (Short 
Papers, Yol. aG- 12, pp 585-588, October 1967. 

E.C. Ba:*be, Linear control systems (Book), 

Interna';! onal TeSt Book Company, Scranton, Pennsylvania, 
pp 262-166, 1963, 

L.D, Hairis, Introduction to feedback system (Book), 

John Wilsy & Sons Inc., New York, pp 161-163, 1961, 

Samuel B^r of sky. Elementary theory of equations (Book), 
The Macl/cillan Company, New York, pp 44-45, 1950, 

J. G, Truxal, Automatic feedback control system 
synthesis (Book), McGraw-Hill Book Company, Inc., 

New York, N.Y. pp 438-440, 1955. 

K. Ogata, State space analysis of control system (Book), 
Prentice Hall, 1967. 

R.D. Gustafson, private correspondence. 


H.M. Power, "a canoniccil matrix for linear discrete- 
time systems, with application to Liapunov function 
generation", submitted to the Journal of the Frajoklin 
Institute. 


P.C, Parks, "Liapunov and the Schur - Cohn stability 
criteria", IEEE Trans, on Automatic Control, Yol, AC-9, 
p 121 ,, January 1964, 



